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Kapitel 1 


EINFUHRUNG 


1.1 DAS BMBF-PROJEKT ILAS 


iLaS steht für intelligentes Laserlicht für kompakte und hochauflösen- 
de adaptive Scheinwerfer. BmBF für Bundesministerium für Bildung 
und Forschung. 


Laut dem Jahresbericht des Deutschen Verkehrssicherheitsrats von 
2011 sind die Unfallzahlen auf deutschen Straßen in den letzten zehn 
Jahren um 25 % gesunken [1]. Gleichzeitig konnten durch die Einfüh- 
rung aktiver und passiver Sicherheitselemente in den Fahrzeugen der 
Insassenschutz und damit auch die Zahl der verunglückten Perso- 
nen im Straßenverkehr deutlich gesenkt werden [1]. Gemessen am 
Fahrzeugaufkommen stellt der überdurchschnittlich hohe Anteil an 
Unfällen bei schlechten Sichtverhältnissen, speziell in der Damme- 
rung und bei Nacht, eine besondere Herausforderung dar [2]. Über 
die Nutzbarmachung der höheren Leuchtdichte und damit verbun- 
denen Auflösungsmöglichkeiten von Laser-Weißlicht ließe sich eine 
noch exaktere Ausleuchtung des Verkehrsraums gewährleisten [2], 


[3]. 


Mit der Lasertechnologie können Leuchtdichten für Weißlicht erzielt 
werden, die voraussichtlich um einen Faktor zwei bis fünf größer 
gegenüber der LED-Technologie von 2014 sind [4]. 
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Mit dem Verbundprojekt sollen die Grundlagen fiir den Aufbau ei- 
nes neuartigen Scheinwerfersystems erforscht werden, das tiber den 
Einsatz eines scannenden Systems im Fern- und Abblendlichtbereich 
eine auf die jeweilige Situation („situativ“) intelligente und anpassbare 
(„variable“) Ausleuchtung des Straßenraums ermöglicht. Zusätzlich 
sollen Licht- bzw. Komfortfunktionen wie etwa Markierungs- oder 
Kennzeichnungsbilder im Nahbereich erarbeitet, im Rahmen von phy- 
siologischen Untersuchungen geprüft und in das Tisch- sowie das 
Fahrmuster zur Validierung im Rahmen von Tests und Probandenstu- 
dien integriert werden. [5] 


Dieses innovative Scheinwerfersystem soll direkt zur Steigerung der 
Verkehrssicherheit beitragen, indem es ein besseres Licht bzw. eine grö- 
Sere wahrgenommene Helligkeit als aktuell auf dem Markt befindliche 
Xenon-Scheinwerfer liefert und gleichzeitig eine höhere Auflösung 
bereitstellen, als das 2013 im Audi A8 in die Serie eingeführte LED- 
System mit einzeln adressierbaren LEDs. Bei erfolgreichem Abschluss 
des Vorhabens und dem Nachweis zur Einsetzbarkeit dieses inno- 
vativen Ansatzes, ist der Grundstein für eine Markteinführung im 
Automobilsektor gelegt. Damit besäße Deutschland gegenüber der 
asiatischen und nordamerikanischen Konkurrenz einen langjährigen 
Wettbewerbsvorteil hinsichtlich des Einsatzes von Laserlichttechno- 
logie für Beleuchtungsfunktionen und Lichtsysteme von Fahrzeugen. 


[5] 


1.2 BLENDFREIES FERNLICHT 


Blendfreies Fernlicht ist eine kameragestützte, dynamische Beleuch- 
tungssteuerungsstrategie, die zielgerichtet Teile der Fernlichtvertei- 
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lung ausblendet, um andere Verkehrsteilnehmer vor Blendung zu 
schiitzen, wahrend dem Fahrer durchgehend die maximale Sichtweite 
geboten wird [6], [7]. 


Die Umgebung der anderen Verkehrsteilnehmer wird ständig mit 
der hohen Lichtintensität des Fernlichts beleuchtet, jedoch ohne die 
Blendung, die normalerweise durch dessen Verwendung verursacht 
würde. Dieses sich ständig ändernde Strahlmuster erfordert komplexe 
Erkennungs-, Verfolgungs- und Steuerungsstrategien, da die Fahrzeu- 
ge, die aus dem Strahl abgeschattet werden müssen, in Bewegung sind. 
Die dynamische Abschattung kann mit beweglichen Schattenmasken 
erreicht werden, die innerhalb des Lichtpfads im Scheinwerfer ver- 
schoben werden. Der Effekt kann auch durch selektives Abdunkeln 
adressierbarer LED-Emitter oder Reflektorelemente erzielt werden, 
eine als Pixellicht bekannte Technik. [6], [7] 


Dieses Prinzip ist in diesem Projekt verbessert worden [3]: statt der 
digital abschaltbaren LED-Emitter ist es nun möglich dem entgegen- 
kommenden Verkehr analog und ohne Sprünge in der Lichtverteilung 
mit einem abgeschatteten Bereich zu folgen, was im folgenden Ab- 
schnitt dargestellt ist. 


1.3 STRUKTUR DER ARBEIT 


In dieser Arbeit werden die thermischen Herausforderungen an ein 
blendfreies Fernlicht untersucht, das mittels mehrerer Laserstrahlen 
und einem Weißlichtkonverter erzeugt wird. Da sich der Aufbau der 
Arbeit am Strahlengang orientiert, ist das optische System des iLaS- 
Scanners im Folgenden dargestellt. 
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1.3.1 DAS OPTISCHE SYSTEM DES ILAS-SCANNERS 


Während dieses Projektes wurden sehr viele optische Konzepte analy- 
siert (siehe dazu [8]). Dies schlägt sich unter anderem in der Vielzahl 
der untersuchten Laserquellen in Kapitel 2 nieder. Schlussendlich hat 
sich jedoch herauskristallisiert, dass ein Scanner mit sechs überein- 
ander angeordneten Spots, der in horizontaler Richtung scannt, der 
sinnvollste Ansatz innerhalb der im Projekt gegebenen Rahmenbe- 
dingungen ist [3], [8]. Dessen Aufbau ist schematisch in Abbildung 
1.1 gezeigt. In diesem Optiksystem fokussiert die Asphäre eine der 
Achsen des Laserstrahls auf den Konverter, wie in Abbildung 1.1(a) 


gezeigt. 


Stellt man sich die Gehäuse der Laserdioden als Kreise vor, die auf 
einer Linie angeordnet sind und aus dem Zentrum der Kreise tritt 
der Laserstrahl aus, ist ersichtlich, dass der minimale Abstand zweier 
Laserstrahlen einem Kreisdurchmesser entspricht. Dies ist für den An- 
wendungsfall als blendfreies Fernlicht nicht eng genug (siehe 1.2(b)). 
Daher müssen die Laserstrahlen gebündelt werden, was in Abbildung 
1.1(b) zu sehen ist. 


Nach der Bündelung wird das Strahlpaket über einen beweglichen 
MOEMS-Spiegel (MOEMS: Micro-Opto-Electro-Mechanical Systems) 
geleitet. Mittels einer weiteren Linse, wie sie schematisch in Abbil- 
dung 1.1(c) als Primäroptik eingezeichnet ist, wird die zweite Achse 
fokussiert. Der Spiegel lässt den Spot-Stapel auf dem Konverter hin 
und her wandern, sodass sich die gewünschte Lichtverteilung, z.B. 
wie in Abbildung 1.2(a), ergibt. 


Der Konverter erzeugt aus dem blauen Laserlicht über Konversion 
weißes Licht. Dieses wird mittels einer Sekundäroptik in das Fernfeld 


projiziert. 
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Die thermisch kritischen Bauteile sind dem Stahlengang folgend: 
1. die Laserdioden selbst, 
2. die Anbindung der Laserdioden an das Scheinwerferkühlsystem, 
3. das Zusammenspiel von Laserdiode und Asphare, 


4. der Konverter. 


oyuro 
-s3unjspungfyens 


(c) MOEMS, Primäroptik, Konverter und Sekundäroptik 


Abbildung 1.1: Das optische System des iLaS-Scanners. a) die Asphäre fokussiert eine 
der Achsen des Lasers auf den Konverter. b) Die Laserstrahlen werden gebündelt, 
d.h. näher zusammen geführt. c) Die gebündelten Strahlen werden über den MOEMS- 
Spiegel abgelenkt und die zweite Achse der Laserstrahlen auf den Konverter mittels 
der als Primäroptik bezeichneten Linse fokussiert. Der MOEMS-Spiegel scannt die 
Laserstrahlen über den Konverter, der aus dem blauen Laserlicht weißes Licht erzeugt. 
Dieses wird über die Sekundäroptik ins Fernfeld abgebildet. Siehe auch [8]. 
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cd] min 


(a) Lichtverteilung mit eng tibereinander angeordeten Laserspots. 


cd] min 


(b) Lichtverteilung mit weiter gespreizten Laserspots. 


Abbildung 1.2: Mégliche Lichtverteilungen zur Veranschaulichung. Weitere Lichtvertei- 
lungen sind in [3] zu finden. 


1.3.2 DIE LASERDIODEN 


Das korrekte thermische Management stellt bei modernen Hochleis- 
tungs-Halbleiterlasern einen Schlüsselfaktor für den reproduzierba- 
ren Betrieb dar. Bei Hochleistungs-Multimodenhalbleiterlasern ist die 
Leistungsdichte höher als bei vergleichsweise schwachen Einzelmo- 
denhalbleiterlasern. Dadurch wird deutlich mehr Verlustwärme pro- 
duziert. Besonders im Dauerbetrieb ist diese Einflussgröße von Bedeu- 
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tung, da durch verschiedene Temperaturen optische sowie geometri- 
sche Anderungen hervorgerufen werden. 


Werden Laserdioden zu heiß, d.h. außerhalb der Spezifikation betrie- 
ben, treten neben schnellerer Alterung, die ebenso bei Betrieb inner- 
halb der Spezifikation auftritt, spontane, unkalkulierbare Ausfälle auf 
[9]. Dies ist für den sicheren Betrieb im Straßenverkehr inakzepta- 
bel. 


Zur Vermeidung dieser Schäden werden in Kapitel 2 zunächst ver- 
schiedene, zur Auswahl stehende Laserdioden analysiert und anhand 
ihres thermischen Widerstands und ihrer Effizienz verglichen. Der 
thermische Widerstand Ry, (auch Wärmewiderstand, Wärmeleitwi- 
derstand) sagt aus, welche Temperaturdifferenz in einem Objekt beim 
Hindurchtreten eines Wärmestromes entsteht [10]. Im Fall der Laserdi- 
ode also, um wie viel höher die Junction-Temperatur im Vergleich zur 
Gehäusetemperatur bei gegebener Verlustleistung ist. Je kleiner der 
thermische Widerstand, desto geringer die Temperaturdifferenz und 


umso kleiner die Junction-Temperatur. 


1.3.3 SCHNITTSTELLE ZUM SCHEINWERFERKÜHLSYSTEM 


Die Gehäusetemperatur und damit, wie eben beschrieben, die Junction- 
Temperatur und damit der sichere Betrieb des iLaS-Scanners, hängt 
maßgeblich von der Anbindung des Lasergehäuses an das Kühlsystem 
ab. Kühlkörper und Lüfter werden von einem assoziierten Partner, 
dem Scheinwerferhersteller ZKW, zur Verfügung gestellt. Daher wird 
in Kapitel 3 die Anbindung des Scanners an die Kühlfläche mittels 
Simulationen untersucht. Dazu wird für jede Laserdiodenvariante 
die maximal zulässige Temperatur ermittelt, die an der Schnittstelle 
zwischen Scanner und Kühlkörper herrschen darf, um den sicheren 
Betrieb noch gewährleisten zu können. Anhand dieser Simulationen 
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wird die Auswahl der fiir das Projekt verwendeten Laserdiode getrof- 
fen. Die Ergebnisse der Simulationen werden mittels Messungen mit 


einer Infrarotkamera am Prototypen verifiziert. 


1.3.4 ZUSAMMENSPIEL VON LASERDIODE UND ASPHARE 


In Abbildung 1.1 ist zu erkennen, dass die Weglänge von Laserdiode 
zu Asphäre wesentlich kürzer ist als der Abstand von Asphäre zum 
Konverter. Damit ergibt sich nach dem Strahlensatz, dass eine kleine 
Änderung der Laserdiodenposition, relativ zur Asphärenposition, gro- 
ße Auswirkungen auf das Lichtbild haben kann. Daher wird in Kapitel 
4 zunächst die thermische Ausdehnung der verwendeten Laserdiode 
simuliert. Im Anschluss wird die Auswirkung der Verschiebung der 
Laserfacette auf die Spotposition berechnet. 


Einen weiteren Einflussfaktor stellt der angelegte Vorwärtsstrom dar. 
Es ist wohlbekannt, dass in Halbleiterlasern die Halbwertsbreite [11] 
und die Strahldivergenz [12] und damit auch die Fokusposition vom 
Vorwärtsstrom abhängt. Daher wird in Abschnitt 4.3 messtechnisch 
untersucht, wie groß der Effekt für das in diesem Projekt verwendete 
optische System im statischen, d.h. mit ausgeschaltetem Scanspiegel, 
sowie im dynamischen Fall, d.h. mit eingeschaltetem Scanspiegel, 
ist. 


1.3.5 DER KONVERTER 


Einer der wichtigsten einschränkenden Faktoren eines Leuchtstoffs 
in einer Anwendung mittels Laser Activated Remote Phosphor, kurz 
LARP, ist der Verlust der Konversionseffizienz bei intensiver Pumpan- 
regung. Die Hauptverluste sind das thermische Quenching der Emis- 


sion, was die Quanteneffizienz bei erhöhten Temperaturen verringert 
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und das nicht-thermische Quenching, das durch die Pumpintensitat 
verursacht wird [13], [4], [14]. Dadurch wird der Anteil an konvertier- 
tem Licht verringert, was neben einer Verringerung des Lichtflusses 
auch eine Farbänderung hervorruft. Dies kann soweit gehen, dass die 
resultierende Farbe nicht mehr die gesetzlichen Vorschriften für die 
Lichtfarbe „weiß“ erfüllt. 


In Kapitel 5 wird die Diffusionsnäherung der Strahlungstransport- 
gleichung zum ersten Mal zur Simulation eines dreidimensionalen, 
nicht rotationssymmetrischen, keramischen Leuchtstoffs eingesetzt. 
Außerdem wird dieser Leuchtstoff erstmalig mittels der Diffusionsnä- 
herung mit einer zeitlich und örtlich veränderlichen Laserlichtquelle 
simulatorisch untersucht. 


1.3.6 DIE KONVERTERMESSUNGEN 


Die Simulationsergebnisse, die mittels der Diffusionsnäherung berech- 
net wurden, werden zunächst mit thermographischen Messungen und 
im Anschluss mit Leuchtdichtemessungen verglichen. 


Die thermischen Messungen werden mit einer Infrarotkamera durch- 
geführt, die Leuchtdichtemessungen mit einer Leuchtdichtekamera. 
Beide Systeme müssen zunächst kalibriert werden. Dazu wird der 
Emissionskoeffizient des Konverters im Infrarotbereich der Thermo- 
graphiekamera bestimmt. Für die Leuchtdichtekamera wird ein Farb- 
vergleich mit einem Spektrometer durchgeführt. Im Anschluss werden 
die jeweiligen Messergebnisse den Simulationsergebnissen gegenüber- 
gestellt. 


Kapitel 2 


THERMISCHE ANALYSE VON 
HALBLEITERLASERN ALS 
PUMPLICHTQUELLE 


Es war H.J. Round, der im Jahre 1907 als erster das Phänomen der 
Lichterzeugung durch elektrische Anregung eines Festkörpers, eines 
SiC-Kristalls, beobachtete, als er an Radiodetektoren forschte [15], 
[16]. Es wurde bereits damals erkannt, dass es sich um kaltes Licht 
handelt, bei dessen Erzeugung sich der Kristall nicht merklich er- 
wärmt. Diese Entdeckung geriet jedoch in Vergessenheit und wurde 
erst im Jahre 1921 wiederentdeckt [16]. Es dauerte weitere 25 Jahre, 
bis der Effekt von Lehovec befriedigend erklärt wurde [16]. Ende der 
Fünfzigerjahre erwartete man noch, dass sich mit LEDs ein Quan- 
tenwirkungsgrad von 100 % erreichen lässt, was bedeutet, dass jedes 
eingebrachte Elektron auch ein Photon erzeugt [16]. Erste Erfolge in 
der LED-Produktion stellten sich mit der Verwendung von GaS- und 
GaP-Kristallen ein. 


Der erste Halbleiterlaser wurde 1962 in verschiedenen Forschungs- 
gruppen nahezu gleichzeitig hergestellt [16]. Diese konnten jedoch nur 
bei kryogenischen Temperaturen betrieben werden [17]. 


Heute ist es zwar möglich Halbleiterlaser bei Raumtemperatur zu 
betreiben, jedoch stellen die Anforderungen in einem Scheinwerfer 
die Dioden vor Herausforderungen. Die Sperrschicht- oder Junction- 
Temperatur (Ty) einer Laserdiode darf nicht über einen kritischen Wert 
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ansteigen, da sonst spontane, zerstörerische Effekte auftreten. Betreibt 
man den Laser unter jener kritischen Temperatur, so tritt lediglich die 
normale Altersdegradierung auf. Diese läuft bei hohen Temperaturen 
schneller ab als bei niedrigen. Die Alterung ist in erster Linie auf die 
Vergrößerung von Fehlstellen im Kristall zurückzuführen [9]. Diese 
Stellen nehmen nicht mehr an der Lichterzeugung teil [9]. Es entstehen 
strahlungslose Übergänge [9]. Somit ist abzuwägen, wie aufwändig 
das Kühlsystem im Verhältnis zur erwarteten Lebensdauer ausfallen 
soll oder muss. 


Bei dieser Betrachtung gilt es auch zu bedenken, dass die heutigen 
Laserdioden (im Jahr 2022) weit von dem in den Fünfzigerjahren 
erwarteten Wirkungsgrad von 100 % entfernt sind [18]. Die blauen 
Laserdioden, die in diesem Projekt Verwendung finden, haben eine 
typische Effizienz von ca. 20-25% [18]. Es ist zu berücksichtigen, 
dass der Wirkungsgrad temperaturabhängig ist [18]. Demnach gilt: 
Je mehr Laserleistung abgerufen werden soll, desto mehr wird sich 
die Sperrschicht erwärmen. Das bedeutet, dass auch hier eine Abwä- 
gung bzgl. der zu erwartenden Lebensdauer und der gewünschten 
Laserleistung und der damit verbundenen Kühlleistung zu treffen 
ist. Diese Betrachtungen werden in den folgenden Abschnitten disku- 
tiert. 


2.1 SIMULATION DES THERMISCHEN WIDERSTANDES 
VERSCHIEDENER HALBLEITERLASER 


Der thermische Widerstand R;n (auch Wärmewiderstand oder Wärme- 
leitwiderstand) sagt aus, welche Temperaturdifferenz in einem Objekt 


12 


SIMULATION DES Ry, VERSCHIEDENER HALBLEITERLASER 


beim Hindurchtreten eines Wärmestromes entsteht [10]. Im Fall der La- 
serdiode also, um wie viel höher die Junction-Temperatur im Vergleich 
zur Gehäusetemperatur bei gegebener Verlustleistung ist. Je kleiner 
der thermische Widerstand, desto geringer die Temperaturdifferenz 
und umso kleiner die Junction-Temperatur und desto besser für die 
Kühlung. 


R = = 
i Pe — Popt Pin 


Tj-T AT K 

J C = | (2.1) 
Dabei sind (P.ı) die elektrische Leistung, (Popt) die optische Leis- 
tung und (Pin) die Wärmeleistung oder auch Verlustleistung, Ty die 
Junction- oder Sperrschichttemperatur und Tc die Case- oder Gehäu- 
setemperatur der Laserdiode. 


Die Anforderungen an die Laserquelle haben sich im Laufe des Pro- 
jekts vielfach geändert (siehe dazu [8]). Schlussendlich hat sich heraus- 
kristallisiert, dass der Laserscanner sechs Kanäle haben soll, welche 
jeweils zwei Watt optische Leistung bei einer Case-Temperatur von 
85 °C liefern können sollen. Die Case-Temperatur (Tc) ist jene Tempe- 
ratur, die am Gehäuse der Laserdiode vorherrscht. 


Bei der Firma OSRAM stehen verschiedene Halbleiterlaser in verschie- 
denen Gehäusen zur Verfügung. Während der Anfertigung dieser 
Arbeit sind die Laserdioden im TO56-Gehäuse, im Multi Die Packa- 
ge, im SMD-Gehäuse und im TO-90 Gehäuse auf ihre Tauglichkeit in 
thermischer Hinsicht untersucht worden. Da in diesem Projekt die 
Grenzen des Machbaren betreffend der Temperaturen ausgelotet wer- 
den sollen, ist bei allen Betrachtungen der thermische Widerstand 
bei der maximal zulässigen Temperatur zu ermitteln und bei allen 
folgenden Berechnungen zu verwenden. Maximal zulässig bedeutet 
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in diesem Fall, die maximale Temperatur, die einen sicheren Betrieb 
über die geforderte Lebensdauer zulässt. 


Da in Datenblättern typischerweise die thermischen Widerstände bei 
Raumtemperatur angegeben werden, muss, gerade wegen der relativ 
hohen Temperaturen, die hier angestrebt werden, für eine korrekte 
Betrachtung auch die Temperaturabhängigkeit der Warmeleitfahig- 
keit der verwendeten Materialien berücksichtigt werden, um eine 
Abschätzung zu erhalten, die möglichst applikationsnah ist. Um dies 
aufzuzeigen wird in Abschnitt 2.1.5 der thermische Widerstand bei 
verschiedenen Gehäusetemperaturen simuliert und das Verhalten gra- 
fisch dargestellt. 


Für die Bestimmung des Rın wurden Simulationen der einzelnen La- 
serdiodenvarianten durchgeführt. Die Randbedingungen der Berech- 
nungen und die Simulationsergebnisse der untersuchten Laserdioden 
werden im folgenden diskutiert. 


2.1.1 SIMULATIONSPARAMETER 


Die Randbedingungen in den Unterabschnitten dieses Kapitels sind 
sehr ähnlich, weshalb sie hier allgemeingültig zusammengefasst dar- 
gestellt werden sollen. Abweichungen und Ergänzungen sind in den 


einzelnen Kapiteln zu finden. 


Da es sich um eine sogenannte interne Simulation handelt, bei der nur 
die Wärmeleitung im Festkörper untersucht wird und ohne Strahlung 
simuliert wird, ist der Berechnungsraum stets exakt so groß wie das in 
den Unterabschnitten dargestellte CAD-Modell. Es wurde die interne 
Studie mit nur in Festkörpern stattfindender Wärmeleitung als Ana- 
lysemethode gewählt, da die Luftzirkulation um die Laserdiode nur 
einen vernachlässigbaren Effekt auf die Gesamtkühlung hat. Gleiches 
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Tabelle 2.1: Allgemeine Randbedingungen der thermischen Simulationen. 


Parameter Einstellung 
Simulationssoftware: SolidWorks 
Flow Simulation 2015 SP5 
Flow Simulation 2017 SP5 


Studienart: Intern 
Warmeleitung in Festkörpern: Ein 
Wärmeleitung nur in Festkörpern: Ein 

Strahlung: Aus 
Zeitabhängig: Aus 
Gravitationseffekte: Aus 

Drehung: Aus 
Strömungsart: nicht anwendbar 
Fluss mit hoher Machzahl: Aus 
Standardrauheit: 0 ym 


Standardzustand der Außenwand: Adiabatische Wand 


gilt für die IR-Strahlung, da die Komponenten aus blankem Metall 
mit einem Emissionskoeffizienten von weniger als 0,1 gefertigt sind. 
Daher sind Strahlung, Gravitationseffekte, Strömung und Fluss mit 
hoher Machzahl ausgeschaltet. 


Das Netz ist so ausgelegt, dass die kleinste Struktur in der Simulati- 
on, der Halbleiter in der Laserdiode, mit mindestens vier Zellen in 
seiner Höhe aufgelöst wird. Ein Beispielnetz ist in Abbildung 2.1 zu 
sehen. 


Die Materialdaten sind, soweit es sich um Standardmaterialien handelt, 
aus der SolidWorks Flow Simulation Materialdatenbank entnommen. 
Die Daten sind in Anhang B zu finden. 
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SIR x 
(a) Gesamte Laserdiode ohne Kappe. (b) Vergrößerung des Chip-Bereichs. 


Abbildung 2.1: Netz für die Simulation der TO90-Laserdiode. 


Als Wärmesenke nach außen ist ein thermischer Kontaktwiderstand 
von 10000 W - m~! - K~? realitätsnah: Dieser Kotaktwiderstand ent- 
spricht einem thermischen Interfacematerial (TIM) mit einer Wärme- 
leitfahigkeit von 2 W - m`! . K7! mit einer Dicke von 200 pm. Außer- 
dem stimmt dieser Wert gut mit Literaturangaben zu metallischen 
Kontakten überein (siehe [19], [20], [21]). Um die Vergleichsgröße der 
Simulationen, die Sperrschichttemperatur der Laserdiode, über die 
verschiedenen Berechnungen möglichst konstant zu halten, wurden 
die Wärmesenken-Temperaturen der untersuchen Varianten entspre- 
chend angepasst. 


Als Konvergenzkriterium ist stets die globale, maximale Festkörper- 
temperatur gewählt. Hierbei kann SolidWorks das Analyseintervall 
und das Abbruchkriterium selbst festlegen oder manuell gegengesteu- 
ert werden. Bei internen Studien mit nur in Festkörpern stattfindender 
Wärmeleitung hat sich dies als nicht nötig herausgestellt. Für exter- 
ne Studien mit Strahlung und Luftstrom kann es nützlich sein, das 
Analyseintervall und das Abbruchkriterium manuell festzulegen, um 
Berechnungszeit zu sparen. 
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Um die Verlustleistung der Laserdioden-Varianten berechnen zu kön- 
nen, muss das Verhalten der Vorwärtsspannung und der emittierten 
Laserleistung über Temperatur und Vorwärtsstrom bekannt sein. Da in 
diesem Projekt ein Stromtreiber zum Einsatz kommt, wird der Strom 


als konstant angenommen. 


Poy = Peı(I,Ty) — Popt(1,Tj) (2.2) 
Pr, = U(LT;) -I— Popt(L Ty) (2.3) 


Die fraglichen Werte können für verschiedene Sperrschichttemperatu- 
ren aus dem Datenblatt der Laserdioden (siehe Anhang A) entnommen 
und mittels eines Tabellenkalkulationsprogramms nutzbar gemacht 
werden. 


Als Beispiel soll hier eine Vergleich der Verlustleistungen bei verschie- 
denen, willkürlichen Betriebsbedingungen der TO90 Laserdiode vor- 
genommen werden. Sie soll zum einen bei Ip = 2,0 A und Tc = 25°C, 
zum anderen bei Ip = 2 0A und Tc = 85°C und zum dritten bei 
Ip = 1,0 A und Tc = 25°C betrieben werden. Man erkennt in Tabelle 
2.2, dass die elektrische Leistung und damit auch die Vorwärtsspan- 
nung sowie die Effizienz 7 der Laserdiode bei hohen Temperaturen 
abnimmt. Ebenso ist ersichtlich, dass die TO90 Laserdiode bei hohen 
Strömen effizienter ist als bei niedrigen. 


Unter den genannten Randbedingungen und der aufgezeigten Me- 
thode zur Berechnung der Verlustleistung können nun thermischen 
Widerstände Laserdioden simuliert werden. Im Anschluss folgt eine 
Zusammenfassung der Simulationsergebnisse und deren Diskussi- 
on. 
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Tabelle 2.2: Vergleich der elektrischen, optischen und thermischen Leistung sowie der 
Effizienz y der TO90 Laserdiode bei verschiedenen Betriebsbedingungen, 
berechnet aus den Datenblattangaben nach Gleichung 2.3. 


IE[A] TePC] PoplW] PalW] PalW] 7 (Popt/Pei) 


2,0 25 2,85 9,16 6,31 31% 
2,0 85 2,16 8,36 6,20 26% 
1,0 85 0,86 3,79 2,93 23% 


2.1.2 LASERDIODE IM TO56 GEHÄUSE 


Die Diode im TO56-Gehäuse ist die erste Laserdiode, die jemals für ein 
LARP-Projekt im Automobilbereich eingesetzt wurde. Es ist folglich 
der logische, erste Schritt, zu versuchen, diese wohlbekannte Laserdi- 
ode auch für iLaS zu verwenden. 


Simulationssetup 


e Wie aus Abbildung 2.2 hervorgeht, sind die Wärmeleitfähig- 
keiten der verwendeten Materialien temperaturabhängig. Das 
bedeutet, dass auch der thermische Widerstand von der Tempe- 
ratur abhängt. Daher sind die hintere und äußere Flächen des 
Gehäuses als Wärmesenke applikationsnah auf 85°C gelegt. 


e Die Diode liefert eine optische Leistung von 1 W. Das bedeutet 
eine thermische Verlustleistung von 3,5 W bei Tc = 85°C. 


e Die epitaktischen Schichten sind in dieser Simulation nicht be- 
rücksichtigt. Da diese recht dünn sind, haben sie nur einen gerin- 
gen Einfluss auf den thermischen Widerstand (siehe Gleichung 


(2.1)). 
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(a) Temperaturabhängige Wärmeleitfähig- (b) Temperaturabhängige Wärmeleitfähig- 
keit A von 6H-SiC [22]. keit A von GaN [23]. 


Abbildung 2.2: Wärmeleitfähigkeiten zur Simulation des Ry, der TO56 Laserdiode bei 
verschiedenen Gehäusetemperaturen (Tc). 


Die Laserdiode im TO-56-Gehäuse leistet optisch bei Raumtemperatur 
maximal zwei Watt. Die Simulationen zeigten, dass die Junctiontem- 
peratur bei Tc = 85°C und Popt = 2W so hoch ist, dass spontane 
Schäden drohen (s. [24]). Um das zu verhindern und die geforder- 
te Lebensdauer zu erreichen, muss Pop: auf ein Watt reduziert wer- 
den. 


Die simulierte Temperaturverteilung in Abbildung 2.3 zeigt, dass die 
Temperatur auch bei nur einem Watt optischer Leistung fast 150 °C 
beträgt. Das bedeutet, dass die Laserdiode, unter zuvor genannten 
Randbedingungen, am äußersten Limit betrieben würde, da 150 °C 


nach Datenblattangaben die maximale Sperrschichttemperatur ist 
[24]. 
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150.00 
140.71 
131.43 
122.14 
112.86 
103.57 
94.29 


85.00 
Temperatur [°C] 


(a) Das Simulationsmodell der TO56 (b) Ergebnis der thermischen Simulation: 
Laserdiode. Ty ~ 150°C 


Abbildung 2.3: Aufbau und Ergebnis der Simulation. 


Der thermische Widerstand betragt laut Simulation und nach Glei- 
chung (2.1) circa 15K - W71. Da die epitaktischen Schichten in der 
Simulation nicht berücksichtigt sind, ist dies eine leichte Unterschät- 


zung. 


Rin(Tc = 85 °C) = 2 = 15 

Dass die Leistung der Laserdiode auf ein Watt beschränkt ist, bedeu- 
tet, dass pro Kanal zwei dieser Dioden Verwendung finden müssten. 
Dies wiederum impliziert, dass zunächst zwei Laserspots auf dem 
Phosphor über den MOEMS-Spiegel genau aufeinander gelegt werden 
müssten, um die Leistung von zwei Watt pro Laserspot zu erreichen. 
Diese Prozedur müsste sechsmal wiederholt werden und diese sechs 
Spots müssten, um ein gutes Ausblendeverhalten zu generieren, alle 
in einer vertikalen Linie angeordnet sein. Da dies nur über sehr kom- 
plexe Optiksysteme möglich ist, wenn überhaupt, wurde diese Lösung 
verworfen. 
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2.1.3 LASERDIODE IM SMD PACKAGE 


Der Chip der SMD-Diode ist eine Neuentwicklung und liefert unter 
normalen Betriebs-Bedingungen bis zu drei Watt optischer Leistung. 
Das Gehäuse basiert auf einem Keramiksubstrat, das mit einem Glas- 
deckel verschlossen und somit der Chip vor Umwelteinflüssen ge- 
schützt wird. Der mittlere Kontakt des SMD-Gehäuses ist elektrisch 
isoliert. Dies bedeutet, dass eine direkte thermische Anbindung an 
den Kühlkörper mehrerer SMD-Laserdioden, beispielsweise durch ein 
Kupferinlay in der Leiterplatte, möglich ist. Daher wird in der Simu- 
lation zur Bestimmung des thermischen Widerstandes angenommen, 
dass die Wärmeabfuhr nur über den mittleren, potenzialfreien Kontakt 
erfolgt. 


Simulationssetup 
e Mittlerer Kontakt als Wärmesenke auf 85°C. 


e Die Diode liefert eine optische Leistung von 2 W. Das bedeutet 
eine thermische Verlustleistung von 5,7 W bei Tc = 85°C. 


e Die Materialdaten für Submount und Die (grün und blau in 
Abbildung 2.5(a)) sind dieselben wie in Abschnitt 2.1.2 in Abbil- 
dung 2.2. Ergänzend zeigt Abbildung 2.4(a) die temperaturab- 
hängige Wärmeleitfähigkeit des Package-Substrats. 


Die epitaktischen Schichten sind in dieser Simulation als dünne 
Schicht berücksichtigt. Das bedeutet, sie liegen im Modell nicht 
als Volumenkörper vor, sondern als thermischer Kontaktwider- 
stand zwischen Submount und Die (grün und blau in Abbildung 
2.5(a)). Der temperaturabhängige thermische Kontaktwiderstand 
ist in Abbildung 2.4(b) dargestellt. 
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(a) Temperaturabhängige Warmeleitfahig- (b) Temperaturabhängiger thermischer 
keit des SMD-Package-Substrats [25]. Kontaktwiederstand der Epitaxie [25]. 


Abbildung 2.4: Wärmeleitfähigkeiten zur Simulation des Ry, der SMD Laserdiode bei 
verschiedenen Gehäusetemperaturen (Tc). 


Das Ergebnis in Abbildung 2.5(b) zeigt, dass die Sperrschichttempera- 
tur Tj fast 160 °C beträgt. Damit würde die Laserdiode, unter zuvor 
genannten Randbedingungen, über dem thermischen Limit betrie- 
ben. Der thermische Widerstand beträgt laut Simulation und nach 
Gleichung (2.1) circa 13K - W-!. Da nur der mittlere, potentialfreie 
Kontakt der Laserdiode in der Simulation zur Entwärmung Verwen- 
dung findet, ist dies eine leichte Überschätzung. 


AT K 
Ra (Te = 85°C) = Pa X By 
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> 150.63 
; 141.25 
131.88 
122.50 
113.13 
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94,38 
85.00 
Temperatur [°C] 
(a) Das vereinfachte Modell der SMD (b) Ergebnis der thermischen Simulation: 
Laserdiode ohne Glasdeckel. Ty ~ 160°C. 


Abbildung 2.5: Aufbau und Ergebnis der Simulation. 


2.1.4 MULTI-DIE-PACKAGE 


Das Multi-Die-Package (MDP) vereint mehrere Halbleiterlaser in ei- 
nem Gehäuse. Die Firma OSRAM bietet ein Gehäuse an, in dem 20 der 
im vorigen Abschnitt im SMD-Gehäuse untersuchten Laserdioden in 
einer (5 x 4)-Matrix angeordnet sind. Insoweit gelten die dort getätig- 
ten Aussagen analog. 


Simulationssetup 
e Unterseite des Kupfergehäuses als Wärmesenke auf 65 °C. 


e Jede Diode liefert eine optische Leistung von 2 W. Das bedeutet 
eine thermische Verlustleistung von 5,7 W bei Tc = 65°C, damit 
Ty < 150°C. 


e Die Materialdaten sind dieselben wie in Abschnitt 2.1.3 


Das Ergebnis in Abbildung 2.6 zeigt, dass die maximale Temperatur 
fast 145°C beträgt. Das bedeutet, dass die Laserchips, unter zuvor 
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(a) Das OSRAM ,,PLPM4 450“ Multi-Die- (b) Ergebnis der thermischen Simulation: 
Package. Ty © 145°C. 


Abbildung 2.6: Aufbau und Ergebnis der Simulation. 


genannten Randbedingungen, am äußersten Limit betrieben würden. 
Der thermische Widerstand für jeden Chip beträgt laut Simulation 
und nach Gleichung (2.1) 

AT K 


Rwy(Tc = 65°C) = — ax 14 —. 
th( C 65 C) Pu W 


Das optische Konzept für das MDP ist ein anderes als für die restli- 
chen Lichtquellen. Es wäre ein Spiegel zum Einsatz gekommen, der 
den Laserstrahl sowohl in vertikaler als auch in horizontaler Richtung 
ablenkt. Der große Nachteil dieses Systems ist, dass solche Spiegel 
zumindest in einer Achse resonant schwingen. Damit ist die Aufent- 
haltswahrscheinlichkeit am Rand des Leuchtstoffs wesentlich höher 
als in der Mitte. Dies widerspricht der Lichtverteilung eines typischen 
automobilen Fernlichts, bei der das Maximum der Intensität in der 
Mitte des Konversionselements benötigt wird. 
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2.1.5 LASERDIODE IM TO90 GEHÄUSE 


Die Firma OSRAM bietet die Laserdiode aus dem SMD Package aus 
Abschnitt 2.1.3 auch im hier untersuchten TO90-Gehäuse an. Daher 
gelten auch hier die im genannten Abschnitt getätigten Aussagen 
analog. 


Simulationssetup 


e Hintere und äußere Flächen des Stahlgehäuses als Warmesenke 
auf 85°C. 


e Die Diode liefert eine optische Leistung von 2 W. Das bedeutet 
eine thermische Verlustleistung von 5,7 W bei Tc = 85°C. 


e Die Materialdaten sind dieselben wie in Abschnitt 2.1.3. 


Die Temperaturverteilung in Abbildung 2.7(b) zeigt, dass die maxi- 
male Temperatur fast 165°C beträgt. Damit würde die Laserdiode, 
unter zuvor genannten Randbedingungen, über dem äußersten Limit 
betrieben. Der thermische Widerstand beträgt laut Simulation und 
nach Gleichung (2.1) 


Die Simulation kann bei verschiedenen Gehäusetemperaturen wieder- 
holt werden. Daraus ergibt sich der temperaturabhängige thermische 
Widerstand der TO90 Laserdiode aus Abbildung 2.8, welcher auch 
in allen folgenden Simulationen Verwendung finden wird. Um ein 
vorhersagefähiges Ergebnis zu erhalten, muss also bei Betriebsbedin- 
gungen simuliert werden, da der Ry, im Datenblatt üblicherweise bei 


Raumtemperatur angegeben wird. 


25 


THERMISCHE ANALYSE VON HALBLEITERLASERN ALS PUMPLICHT 


165.00 
155.00 
145.00 
135,00 
125.00 
115.00 
105.00 
95,00 


85.00 
Temperatur [°C] 


(a) Das Simulationsmodell der TO90 (b) Ergebnis der thermischen Simu- 
Laserdiode. lation: Tj © 165°C. 


Abbildung 2.7: Aufbau und Ergebnis der Simulation. 
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Abbildung 2.8: Temperaturabhangiger thermischer Widerstand der TO90 Laserdiode. 
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2.1.6 ZUSAMMENFASSUNG 


Die Simulationen der zur Verfügung stehenden Laserdioden haben 
ergeben, dass die thermischen Widerstände der einzelnen Varianten 
nicht wesentlich variieren, wie in Tabelle 2.3 zusammengefasst. Ledig- 
lich die Diode im TO56-Gehäuse sticht heraus, da sie bei Erzeugung 
von 1W optischer Leistung 3,5 W thermische Verlustleistung erzeugt. 
Die drei anderen Varianten generieren bei 2W Ausgangsleistung nur 
5,7W Abwärme. Betrachtet man die verbleibenden drei Laserdioden 
einzeln, ohne Halter und Anbindung an ein Kühlsystem, wie in die- 
sem Kapitel geschehen, könnte man aus thermischer Sicht jede der 
Varianten wählen. 


Tabelle 2.3: Thermische Widerstände, optische Leistungen, Verlustleistungen und Effi- 
zienzen der verschiedenen Laserdioden. 


Laser Rm[K-W1] Pop[W] PmI[W] 17 = Popt/Pa 


TO56 15 1,0 3,5 0,22 
SMD 13 2,0 5,7 0,26 
MDP 14 2,0 5,7 0,26 
TO90 14 2,0 5,7 0,26 


Da der thermische Widerstand temperaturabhängig ist, muss bei Be- 
triebsbedingungen simuliert werden, um ein vorhersagefähiges Ergeb- 
nis zu erhalten, vor allem, da der Ry, im Datenblatt üblicherweise bei 
Raumtemperatur angegeben wird. 


Ob die Laserdioden auch in Verbindung mit Halter und Kühlsystem 
untereinander vergleichbare thermische Eingenschaften haben, wird 
im kommenden Kapitel untersucht. 
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Kapitel 3 


SCHNITTSTELLE ZUM 
SCHEINWERFERKUHLSYSTEM 


3.1 SIMULATION DER SCHNITTSTELLEN FUR DIE 
VERSCHIEDENEN LASERDIODEN 


Dieses Kapitel behandelt die Schnittstelle des von OSRAM entwickel- 
ten Laser-Scanner-Moduls zu dem von ZKW beigesteuerten Kühlsys- 
tem, das das Einhalten der im vorigen Kapitel berechneten maximalen 
Junction-Temperaturen garantieren soll. ZKW ist u. a. Scheinwerfer- 
hersteller und war assoziierter Partner im iLaS-Projekt. Zu einer ge- 
lungenen Schnittstellenbeschreibung gehört zum einen die über die 
Schnittstelle abzuführende Leistung sowie die Leitungsverteilung auf 
die einzelnen Verbraucher und zum anderen die maximal zulässige 
Temperatur der Schnittstelle. Außerdem sollte ein CAD-Modell zur 
Verfügung gestellt werden. 


Um die Temperatur der Schnittstelle ermitteln zu können, wurden 
für die Laserdiodenoptionen thermische Simulationen durchge- 
führt. 


Die Aussagen zum Simulationsaufbau aus Kapitel 2.1.1 gelten ana- 
log. 
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3.1.1 SCHNITTSTELLE FUR DIE LASERDIODE IM TO56-GEHÄUSE 


Wie in Abschnitt 2.1.2 beschrieben, sind fiir die geforderten 2W 
pro Kanal zwei TO56 Laserdioden nötig. Die Anordnung der Di- 
oden und das Ergebnis der Simulation ist in Abbildung 3.1 darge- 
stellt. 


SIMULATIONSSETUP 


e Modulrückseite als Wärmesenke auf 85 °C. 


e Jede Diode liefert eine optische Leistung von 1W. 


SIMULATIONSERGEBNISSE 


° To ~ 115°C 
e = Tj = 150°C (siehe 2.1.2) 
e = Tint < 70°C 


Um die Sperrschichttemperatur (Ty) von unter 135°C zu halten, darf 
die Modulrückseite eine Temperatur von 70°C nicht überschreiten. 
Die Junction-Temperatur von 135°C darf nicht überschritten werden, 
um die in diesem Projekt geforderten Lebensdaueranforderungen 


garantieren zu können. 
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(a) Anordnung der zwölf TO56 Laserdi- 
oden (b) Ergebnis der thermischen Simulation 


Abbildung 3.1: Thermische Simulation von zwölf TO56 Laserdioden. Geometrische 
Anordnung zur Kombination je zweier Laserdioden in einem Kanal. 


3.1.2 SCHNITTSTELLE FÜR DIE LASERDIODE IM MDP 


Um auch für diese Lichtquelle den optimalen thermischen Operations- 
punkt zu finden, wurde sie simuliert. Da das MDP auf einer Leiter- 
platte moniert werden kann, wurde diese in verschiedenen Varianten 
in den Simulationen berücksichtigt (siehe Abbildung 3.2). Zum einen 
als Metallkernplatine, zum anderen als FR-4-Platine mit einem Loch, 
sodass die Rückseite des MDP direkt auf dem Kühlkörper anliegt, die 
Kontakte jedoch mittels der Leiterbahnen auf der Platine kontaktiert 
werden können. Da die Treiberbausteine, die den Strom für das MDP 
liefern, auf Grund von EMV-Überlegungen, möglichst nahe am MDP 
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sitzen sollen und diese eine nicht zu vernachlassigende Verlustleistung 
erzeugen, wurden auch sie in der Simulation beriicksichtigt (siehe Si- 
mulationssetup). Sie sind im Fall der FR-4-Platine mittels Kupferinlays 
(orange-farbige Zylinder in Abbildung 3.2(b)) mit dem Kühlkörper 


verbunden. 


(a) MDP auf Metallkernplatine. (b) MDP eingebettet in FR-4-Platine. 


Abbildung 3.2: Simulationsmodelle zur Analyse der verschiedenen Montagemöglich- 
keiten des MDP. 


SIMULATIONSSETUP 


e Rückseite Kupfer-Kühlkörper: 65 °C 


e TIM (1,5W-m-!-K~!, 144pm) zwischen Kühlkörper (blau) 
und MCPCB 


e Dielektrikum (3W-m~!-K~!, 381m) zwischen Kupferkern 
(1mm, grau) und Top-Layer (70 um, Orange) 


e 4 Treiberbausteine mit je 4W thermischer Leistung 
e MDP mit 20 x 5,7 W thermischer Leistung 
e Gehäuse und Innenleben des MDP aus Kupfer 
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e Die Materialdaten sind dieselben, wie in Abschnitt 2.1.4, wenn 
hier nicht anders angegeben. 


(a) MCPCB: Schnitt durch (b) FR-4: Schnitt durch die (c) FR-4: Schnitt durch die 


‚85.00 
183.46 
8192 
8038 
7885 
7731 
75.77 
74.23 
7269 
7115 
69.62 
63.06 
6654 
65.00 


? 
F 


die obere Kupferla- obere Kühlkörperebene obere Platinenebene 
ge zur Analyse der zur Analyse der MDP- zur Analyse der FET- 
MDP und FET Gehäuse- Gehäuse-Temperatur. Gehäuse-Temperatur. 
Temperatur. 


Abbildung 3.3: Die Ergebnisse zeigen einen deutlichen Vorteil der FR-4-Platine gegen- 
über der Metallkernplatine. 


SIMULATIONSERGEBNISSE 


Die Ergebnisse in Abbildung 3.3 zeigen einen deutlichen Vorteil der 
FR-4-Platine gegenüber der Metallkernplatine, da das MDP direkt 
an den Kühlkörper angebunden werden kann. Wie in Kapitel 2.1.4 
beschrieben, benötigt man den Großteil des Lichtes in der Mitte des 
Konverters. Am Rand kann demnach in diesem speziellen Optiksys- 
tem die Laserleistung reduziert werden. Nimmt man eine über ei- 
ne Periodendauer gemittelte Leistung von rund 60% (Py, ~ 78W) 
an, so erhält man mit dem in Kapitel 2.1.4 berechneten thermischen 
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Widerstand eine maximale Kühlkörpertemperatur von 64 °C für die 
Metallkernplatine und 68 °C für die direkte Anbindung an den Kühl- 
körper. 


3.1.3 SCHNITTSTELLE FÜR DIE LASERDIODE IM SMD-GEHÄUSE 


Wie in den Abschnitten zuvor soll auch für das SMD-Package die 
maximale Interface-Temperatur zum Kühlkörper berechnet werden. 
Da die SMD-Laserdiode auf eine Platine gelötet werden muss, sind 
zwei Setups untersucht worden: 


1. Konventionelle Metallkernplatine (siehe Abbildung 3.4(b)), 


2. Metallkernplatine mit Pedestal-Technologie (siehe Abbildung 
3.4(c)). 


Im Vergleich zu Option 1 hat Variante 2 den Vorteil, dass das Dielek- 
trikum an der Stelle des thermischen Pads des Laserdioden-Packages 
ausgefräst und die entstandene Lücke im Galvanikprozess mit Kup- 
fer aufgefüllt wird (siehe Abbildung 3.4). Somit entsteht ein direkter, 
metallischer Pfad vom thermischen Pad des Laserdioden-Packages 
zum Kühlkörper. Nachteil dieser Lösung ist der vergleichsweise hohe 
Preis. 


SIMULATIONSSETUP 


e Aufbau nach Abbildung 3.4(a) 


e Aluminiumkühlkörper als Wärmesenke mit Rückseitentemp. 
85°C 


e Wärmeleitpaste zw. Kühlkörper und Metallkernplatine (MCPCB) 
e MCPCB, Dielektrikum: A = 3W-m7!-K~1!,d = 38 um 
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Kupferenddicke von 75 pm 
2 W optische Leistung pro Laserdiode 


75% duty cycle: (5,2 -0,75)W = 3,9W Verlustleistung pro Laser- 
diode (Reduzierung der Laserleistung am Konverterrand) 


Wärmeleitung nur im Festkörper (keine Konvektion / Strahlung) 


Die Materialdaten sind dieselben wie in Abschnitt 2.1.3, sofern 
hier nicht anders angegeben. 
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< 


(a) Aufbau der Simulation. Schnitt entlang A. 


Substrat 


Anode 1 Therm. F Kathode 


Dielektrikum 


(b) Vereinfachte Schnittdarstellung: Konventionelle 
Metallkernplatine mit durchgehendem Dielek- 
trikum. 


Substrat 


Anode 1 Therm. F Kathode 


Dielektrikum Bl 


(c) Vereinfachte Schnittdarstellung, Pedestal Tech- 
nologie: Dielektrikum unter thermischem Pad 
ausgefräst und aufgekupfert. 


Abbildung 3.4: Unterschiede der konventionellen MCPCB- im Vergleich zur Pedestal- 


Technologie. 
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SIMULATIONSERGEBNISSE 


Die Resultate der thermischen Simulation in Abbildung 3.5 zeigen 
einen Vorteil der Pedestal-Technologie gegenüber der konventionellen 
Metallkernplatine. Während sich auf Grund des höheren thermischen 
Gesamtwiderstands für die klassische Metallkernplatine eine maxi- 
male Kühlkörper-Temperatur von 70°C ergibt, liegt diese bei Verwen- 
dung der Pedestal-Technologie ca. 8°C höher. Die maximale Junction- 
Temperatur dieser Laserdiode liegt bei 135 °C. 


124 

117.50 
111.00 
104.50 
98.00 

91.50 

85.00 

Temperature (Solid) [°C] Temperature (Solid) [°C] 


(a) MCPCB: Die Junction-Temperatur des (b) Pedestal: Die Junction-Temperatur des 
Lasers liegt bei ca 150°C. Lasers liegt bei ca 142°C. 


Abbildung 3.5: Die Ergebnisse zeigen einen deutlichen Vorteil der Pedestal-Technologie 
gegentiber der konventionellen Metallkernplatine. 


3.1.4 SCHNITTSTELLE FUR DIE LASERDIODE IM TO90-GEHAUSE 


Abschließend wird die Interface-Temperatur zum Scheinwerferkühl- 
system für die Laserdiode im TO90 Gehäuse untersucht. Der Aufbau 
der Simulation ist in Abbildung 3.6 gezeigt. 
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SIMULATIONSSETUP 


e Aufbau nach Abbildung 3.6 


e Aluminiumkühlkörper als Wärmesenke mit Rückseitentempera- 
tur 80 °C 


2 W optische Leistung pro Laserdiode 


75% duty cycle: (5,7 -0,75)W = 4,3W Verlustleistung pro Laser- 
diode (Reduzierung der Laserleistung am Konverterrand) 


Die Materialdaten sind dieselben, wie in Abschnitt 2.1.5. 


Abbildung 3.6: Aufbau der Simulation des TO90-Moduls inkl. Linsen. 
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(a) Gesamtmodul mit ausgeblendeten Linsen (b) Vergrößerte Ansicht der Laser- 
und Kappen für die Laserkanäle 4-6. diode des Kanals Nr. 4. 


Abbildung 3.7: Ergebnisse der Thermischen Simulation des TO90-Moduls. 


SIMULATIONSERGEBNISSE 


Abbildung 3.7 zeigt eine Junction-Temperatur von 136°C. Dies ist 
1°C höher, als das Datenblatt der Laserdiode zulässt (siehe [18]). 
Damit ist eine maximale Kühlkörpertemperatur von 79°C einzuhal- 
ten. 


3.1.5 ZUSAMMENFASSUNG 


Die Ergebnisse der in Kapitel 3.1 erstellten Simulationen sind in Ta- 
belle 3.1 zusammengefasst. Je höher die mögliche maximale Interface- 
Temperatur ist, desto kleiner und günstiger kann das angeschlossene 
Kühlsystem ausfallen. Man kann erkennen, dass die Lösung mit SMD- 
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Laserdiode und TO90-Laserdiode gleichwertig sind und die Losungen 
mit MDP und TO56 Lasern dagegen abfallen. 


Wahrend diese Arbeit erstellt wurde, war nicht absehbar, ob es 
einen Laser im SMD-Package geben wiirde. Da fiir den Laser im 
TO90-Gehäuse Prototypen verfügbar waren, wurde diese Variante 
gewählt. 


Tabelle 3.1: Ergebnisse der Interface-Temperatur-Simulationen mit der jeweiligen opti- 
schen Leistung und der Verlustleistung. Je höher die mögliche maximale 
Interface-Temperatur ist, desto kleiner und günstiger kann das angeschlos- 
sene Kühlsystem ausfallen. 


Laser T max, interface [°C] P opt [W] P th [W] 


TO56 70 12 x 1,0 42 
MDP 68 20x 1,2 78 
SMD 78 6x15 26 
TO90 79 6x15 26 


Das Interface zum Kühlsystem kann nun wie folgt beschrieben werden: 
Jede TO90 Laserdiode tragt eine thermische Last von 4,3 W. Am Inter- 
face zwischen Kühlkörper und Lasermodul dürfen bei dieser Last und 
Lastverteilung maximal 79 °C anliegen. Die zur Wärmeübertragung 
nutzbare Fläche beträgt 34,2 cm?. 


3.2 MESSUNG DES TO90-LASERMODULS 


Um die Simulationen aus den vorherigen Kapiteln zu bestätigen, wur- 
de am Prototypen eine Messung mit der Infrarotkamera durchge- 
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führt. In diesem Abschnitt wird zunächst der Versuchsaufbau und 
dessen Durchführung beschrieben. Danach muss die Simulation aus 
dem vorigen Abschnitt unter Messbedingungen wiederholt werden. 
Zuletzt werden die Ergebnisse von Simulation und Messung vergli- 
chen. 


3.2.1 VERSUCHSAUFBAU UND -DURCHFÜHRUNG 


Der Messaufbau ist schematisch in Abbildung 3.8 gezeigt. Um die Tem- 
peratur des TO90-Lasermoduls mit der Infrarotkamera messen zu kön- 
nen, wurde es mit einer Farbe mit einem Emissionskoeffizienten von 
e = 0,95 lackiert, wie in Abbildung 3.9 zu sehen ist. 


Abbildung 3.8: Schematischer Messaufbau mit der IR-Kamera. 
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Abbildung 3.9: Zur Messung mit der IR-Kamera lackiertes TO90-Lasermodul. 


Die Treiberplatine musste zur Stromversorgung genutzt werden. Sie 
bietet keine Möglichkeit zur Messung von Vorwärtsspannung und 
-strom. Daher können die Verluste nur abgeschätzt und eine Konsis- 
tenzprüfung mittels Eingangsstrom- und -spannungs-Messung durch- 
geführt werden. 


Der Treiber stellt den Arbeitspunkt der Laserdioden auf 2,0 A ein. In 
den Umkehrpunkten des Spiegels wird die Stromzufuhr für je 6% der 
Zykluszeit unterbrochen, um Überhöhungen in der Leuchtdichte auf 
dem Konverter zu vermeiden. Der duty cycle beträgt also 88 %. Bei 
einer Case Temperatur von 70°C wird nach den Graphen in Anhang A 
eine Vorwärtsspannung von 4,27 V erwartet. Mit einer erwarteten op- 
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tischen Laserleistung von 2,34 W liegt die prognostizierte Verlustleis- 
tung bei 5,17 W pro Laser. Mit einer Gesamttreibereffizienz von rund 
70% liegt die erwartete Eingangsleistung bei 45 W. 


Mittig in der Adapterplatte ist ein Pt-100 Temperaturfühler einge- 
bracht. Er wird zur Regelung der Peltierelemente tiber den Peltier- 
Controller genutzt. Eine Wasserkühlung sorgt dafür, dass die Peltier- 
elemente nicht überhitzen. Die Adapterplattentemperatur wird auf 
50°C eingestellt. 


Der Steuer-PC gibt mittels Vektor-Box und CANoe Steuersignale 
an die Treiberplatine und per separater Software an den Peltier- 
Controller. 


3.2.2 SIMULATION UNTER MESSBEDINGUNGEN 


Um die Messung mit der Simulation vergleichen zu können, musste 
die Simulation an die geänderten Randbedingungen angepasst wer- 
den. Die Berechnung unterscheidet sich zum vorherigen Kapitel durch 
den geänderten Kontaktwiderstand nach außen und die hinzugekom- 
mene Strahlungskomponente bei Tamb = 25°C, da das Modul, um es 
messen zu können, mit Emissionslack (e = 0,95) lackiert wurde (siehe 
Abbildung 3.9). Es musste zur Messung ein Adapter zwischen Modul 
und Wasserkühlung (siehe Abbildung 3.8) eingesetzt werden und das 
Modul ohne TIM mit diesem Adapter verbunden werden. Der Kontakt- 
widerstand pro Fläche wurde auf 1000 W - m~? - K7! festgelegt, die 
Verlustleistung auf die zuvor abgeschätzten 5,17 W. 


Das Ergebnis der Simulation ist in Abbildung 3.10 gezeigt. Im Ver- 
gleich zur Simulation ohne Strahlung ist ein kleiner Ry, errechen- 
bar. 
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Abbildung 3.10: Ergebnisse der zur Messung passenden thermischen Simulation des 
TO90-Moduls. Die Linsen sind in dieser Simulation kalt, da sie keinen Kontakt zu den 
Haltern aufweisen. 


3.2.3 VERGLEICH VON SIMULATION UND MESSUNG 


In der Simulation und bei der Messung wurde auf dem oberen und 
unteren Balken je ein Oberflächenziel definiert, das die maximale 
Temperatur der gewählten Oberfläche ausgibt. Ein Thermographiebild 
ist in Abbildung 3.11 gezeigt. Tabelle 3.2 vergleicht die Messergebnisse 
mit der Simulation. 


Dass Simulation und Messung an den Balken bis auf ein Zehntel- 
grad übereinstimmen, muss als Zufall gewertet werden. Die Messun- 
sicherheit der Kamera beträgt 2 °C. Die Messunsicherheit des Peltier- 
Controllers mit 0,05 °C ist ebenso vernachlässigbar wie die Schwan- 
kung der Temperatur während der Messung von 0,2°C. Außerdem 
wurde die Verlustleistung der Laserdioden anhand der Datenblatt- 


44 


MESSUNG DES TO90-LASERMODULS 


werte abgeschätzt und nicht gemessen. Die Schätzung der Eingangs- 
leistung der Treiberplatine von 45 W und die Konsistenzmessung von 
(46,10 + 0,24) W stimmen nicht innerhalb des Fehlerintervalls überein. 
Die Abweichung beträgt rund 3 %. 


Tabelle 3.2: Vergleich von Simulation und Messung. Sehr gute Übereinstimmung für 
die „Balken“. Gute Übereinstimmung für die Bodenfläche. 


Einheit Simulation Messung 


max. Temperatur Kanäle 4-6 [°C] 62,69 63 +2 
max. Temperatur Kanäle 1-3 [°C] 62,68 63 +2 
max. Temperatur untere Flache [°C] 59,77 57 +2 


Die Simulation und die Messung weisen leicht unterschiedliche Ma- 
terialstande auf, die keinen signifikanten Einfluss auf die Ergebnisse 
dieses Abschnitts haben. 


Nachdem die thermischen Messungen des TO90-Lasermoduls die 
Simulationen bestätigt haben, werden im nächsten Kapitel zu- 
nächst die Auswirkungen der berechneten Erwärmung der La- 
serdiode auf ihre Geometrie mittels mechanischer Simulationen 
untersucht. 
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BE m 


(a) Fokus auf die „Balken“ . 


(b) Fokus auf den „Boden“ . 


Abbildung 3.11: Messung der Temperaturen des TO90-Moduls. Gut erkennbar sind die 
blauen Rechtecke mit denen der Emissionsgrad auf 95 % festgelegt wird. Die gelben 
Rechtecke sind die Messbereiche. Angeheftet an diese ist die maximale Temperatur der 
jeweiligen Messbereiche. 


46 


Kapitel 4 


ZUSAMMENSPIEL VON LASERDIODE 
UND ASPHARE 


Führt man sich das optische System des iLaS-Scanners vor Augen, 
erkennt man, dass die Weglänge von Laserdiode zu Asphäre um ein 
Vielfaches kürzer ist als der Abstand von Asphäre zum Konverter. 
Damit ergibt sich nach dem Strahlensatz, dass eine kleine Änderung 
der Laserdiodenposition große Auswirkungen auf das Lichtbild haben 
kann. Dies gilt auch für den späteren Zusammenbau von Laserdiode 
und Asphäre in der Fertigung. 


Am einfachsten ist eine Justage in der Fertigung zu bewerkstel- 
ligen, wenn die Laserdiode gerade eben über der Laserschwelle 
und bei Raumtemperatur betrieben wird. Ob dies jedoch den ho- 
hen Ansprüchen an die Genauigkeit der relativen Position von 
Laserdiode und Asphäre entspricht, soll im folgenden untersucht 
werden. 


Dazu behandelt dieses Kapitel simulatorisch die thermische Aus- 
dehnung der TO90 Laserdiode, deren Einflüsse auf die Laserdioden- 
Asphären Unterbaugruppe (Submount) sowie messtechnisch die Ver- 
schiebung der Fokusposition über den Vorwärtsstrom. 
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4.1 SIMULATION DER THERMISCHEN AUSDEHNUNG 
DER TO90 LASERDIODE 


Dass sich Halbleiterlaser wahrend des Betriebs und bei hohen Tem- 
peraturen auf Grund der verschiedenen verwendeten Materialien in 
den drei Raumrichtungen verschieden stark ausdehnen, ist wohlbe- 
kannt (z. B. [26] oder [27]). Zumeist wird versucht, diesem Effekt durch 
die Einbringung verschiedener Materialkombinationen im Gehäuse 
entgegen zu wirken (siehe auch [26] oder [27]). Doch eine leichte Ver- 


formung lässt sich nicht vermeiden. 


Um heraus zu finden, wie sich die in diesem Projekt verwendete La- 
serdiode im TO90-Gehäuse (zur Auswahl siehe Kapitel 3.1.5) verhält, 
muss die Verformung dieses Bauteils simulatorisch untersucht wer- 
den. Die zugehörigen Materialdaten sind in den vorigen Abschnitten 
genannt worden oder finden sich in Tabelle B.1. Es gilt zu eruieren, 
wie sich die Geometrie der Laserdiode im Betrieb bei hohen Gehäu- 
setemperaturen im Vergleich zum Betrieb bei Raumtemperatur än- 
dert. 


Das Simulationssetup ist dasselbe wie in Abschnitt 2.1.5. Für beide in 
Abbildung 4.1 gezeigten Fälle ist die Referenztemperatur 25°C. Das 
bedeutet, dass die Verformung auf Grund der Erwärmung durch den 
Betrieb bei 25°C und durch den Betrieb bei 85°C relativ zur nicht 
betriebenen Diode, die sich auf einer Temperatur von 25°C befindet, 
ermittelt wird. Als thermische Last werden die Simulationsergebnisse 
aus Kapitel 2.1.5 verwendet. 


Wie Abbildung 4.1(a) zeigt, ist allein auf Grund des Betriebs bei 25 °C 
eine Verformung von bis zu 0,6 ym zu erwarten. Wird die Gehäuse- 
temperatur auf 85°C erhöht, was annähernd dem späteren Betriebs- 
punkt entspricht, so ist festzustellen, dass die maximale Ausdehnung 
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ca. 8,5 pm beträgt. Außerdem stellt man eine Verkippung der Facette 
fest. Die Werte sind in Tabelle 4.1 dargestellt. Dies hat zur Folge, dass 
die Emitterposition tiber die Temperatur variiert. 


URES (micron) 
8.500 


1.792 
7.083 
6.375 
5.66? 
4958 
4250 
2542 
2833 
2.228 
za 
0.708 
0,000e+ 000 


(a) Maximale thermische Ausdehnung bei (b) Maximale thermische Ausdehnung 
Tc = 25°C ist ca. 0,6 pm bei Tc = 85°C ist ca. 8,5 um 


URES (micron) 
0.605 
0,555 
0.504 
0.454 

0.404 

0.353 


soe S 


0,303 
0.252 
0.202 
0.352 
omr 
6.050 


Abbildung 4.1: Die Ergebnisse zeigen, dass sich eine deutliche Ausdehnung auf Grund 
des Betriebs und der Betriebstemperatur ergibt. Die Verformung ist in den Abbildungen 
stark tibertrieben dargestellt. 


Fur sich genommen hat dies keine Auswirkung auf die Laserdiode. 
Im Zusammenspiel mit der Asphäre jedoch ergibt sich ein anderes 
Bild: Der Abstand zwischen Asphäre und Emitter bei Tc = 25°C 
ist nicht derselbe wie bei Tc = 85°C (siehe Tabelle 4.1). Die Aus- 
wirkungen auf das Lichtbild sind im folgenden Abschnitt darge- 
stellt. 
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Tabelle 4.1: Ergebnisse der mechanischen Simulation des Verhaltens der TO90 Laserdi- 
ode auf Grund thermischer Ausdehnung. 


Verschiebung der Facette in Strahlrichtung 8,61 um 


Tip der Facette an deren rechter Kante 0,08 ° 
Tip der Facette an deren linker Kante 0,05° 
Tilt der Facette 0,06 ° 


4.2 THERMISCHE AUSDEHNUNG IM TO90 
SUBMOUNT 


Dieser Abschnitt behandelt die Auswirkungen auf das optische Sys- 
tem der im vorigen Abschnitt berechneten Verschiebung und Verdre- 
hung der Laser-Facette. Ahnliche Probleme treten beispielsweise in 
der Kommunikationstechnik bei der Einkopplung in Glasfaserkabel 
[28] oder in Laserdruckern [29] auf. Dabei wird versucht, mit ver- 
schiedenen Komponenten-Materialien [29] oder mit unterschiedlichen 
Lotlegierungen, die Linse und Linsenhalter verbinden [28], die negati- 
ven Auswirkungen der unterschiedlichen Ausdehnungskoeffizienten 


im System zu kompensieren. 


Um die Auswirkungen der unterschiedlichen Ausdehnungskoeffizi- 
enten der in diesem Projekt verwendeten Materialien zu verstehen, 
soll zunächst die Asphare in ihrer nominellen Position verbleiben und 
lediglich die Auswirkungen der Erwärmung der Laserdiode simula- 
torisch untersucht werden. In einem zweiten Schritt wird auch die 
Auswirkung der Verschiebung der Linse auf Grund der Ausdehnung 
ihres Halters simuliert. 
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4.2.1 OPTISCHE SIMULATION - LINSE IN NOMINALPOSITION 


Die optischen Simulationen wurden von Dr. Frank Buchmann (OS- 
RAM GmbH) im nicht sequentiellen Modus des Simulationspro- 
gramms OpticalStudio® durchgeführt. Basierend auf der Strahlendatei 
der Laserdiode beschreibt das Modell den Submount aus Laserdiode 
und asphärischer Linse in Nominalposition und enthält einen De- 
tektorbereich in der Bildebene. Die Projektionslängen der einzelnen 
Optik-Baugruppen liegen, je nach Kanal, im Bereich von 192,3 mm 
bis 225,5 mm. Der Aufbau der Simulation ist exklusive der Detektor- 
fläche, die sich im jeweiligen Bildpunkt befindet, in Abbildung 4.2 
gezeigt. 


Abbildung 4.2: Aufbau der optischen Simulation des TO90-Moduls inkl. Linse in Nomi- 
nalposition. 


Die Aufgabe der Asphäre ist das Fokussieren des Laserspots in Bezug 
auf die fast axis in der Bildebene, d.h. auf den Konverter. Zur Berech- 
nung der Auswirkungen der Erwärmung der Laserdiode wird der 
Quellpunkt in der Simulation nach Tabelle 4.1 geneigt und verscho- 
ben. 
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Auf der Detektorflache der Simulation beobachtet man eine Deplat- 
zierung und eine Defokussierung der projizierten Laserfacette in der 
Bildebene. Abbildung 4.3 verdeutlicht diesen Effekt. Die linke Seite 
zeigt die ideale Position, d.h. ohne thermische Ausdehnung und die 
rechte Seite das Laserprofil unter realen Bedingungen, mit thermischer 
Ausdehnung. Die Detektorfläche für die optische Simulation in der 
Bildebene beträgt 1mm x 1 mm und die Leistungsdichte in W - mm? 


dient als Indikator. 


1,0 Incoherent 18 Incoherent 
Irradiance Irradiance 
7,19 6,11 
6,47 5,50 
5,75 4,89 
5,03 4,28 
o 4,31 o 3,67 
3,59 4 3,06 
2,87 2,45 
2,16 1,83 
1,44 1,22 
0,72 0,61 
0,00 0,00 
-1,0 -1,0 
-1,0 0 1,0 -1,0 o 1,0 


X coordinate value X coordinate value 


Y coordinate value 
Y coordinate value 


Abbildung 4.3: Links: Ideale Position des Bildes. Rechts: Verschiebung und Defokussie- 
rung des Bildes auf Grund der Erwärmung der Laserdiode. 


Axiale Schnitte in der fast und slow axis des Laserprofils veranschauli- 
chen den Einfluss der Temperatur detaillierter (siehe Abbildung 4.4). 
Die effektive Verschiebung des Laserprofils in der fast axis beträgt 
84 um und in der slow axis 98 um. Die Defokussierung auf Grund der 
in axialer Richtung verschobenen asphärischen Linse führt zu einer 
ungewünschten Verringerung der Leistungsdichte. 


Dies bedeutet, dass in diesem optischen System, wenn alle sechs La- 
serdioden betrieben werden, eine Uberlappung und / oder eine Lücke 
zwischen den Laserkanälen am Konverter auftreten würde. Das Ergeb- 
nis wäre ein ,,Zebramuster”, das den Fahrer stören würde, womit die 


Vorteile des Systems nicht mehr gegeben wären. 
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Abbildung 4.4: Links: Axialer Schnitt in der fast axis des Laserprofils (vgl. Abbildung 
4.3) in der Bildebene. Rechts: Axialer Schnitt in der slow axis des Laserprofils (vgl. 
Abbildung 4.3) in der Bildebene. Blau: ursprüngliche Position der Facette; Orange: 
Verschobene Position aufgrund von Erwärmung. 


4.2.2 SIMULATION DER THERMISCHEN AUSDEHNUNG DES TO90 
SUBMOUNTS 


So, wie sich die Laserdiode selbst auf Grund ihrer Eigenerwärmung im 
Betrieb verformt, ist auch der Linsenhalter diesem Effekt unterworfen. 
Der Simulationsaufbau der Baugruppe ist in Abbildung 4.5 gezeigt. 
Als Last dient die thermische Simulation aus Kapitel 2.1.5, erweitert 
um die in Stahl gefasste Linse und den stählernen Linsenhalter. Auch 
die Materialdaten sind dieselben. 


Abbildung 4.6 zeigt das Ergebnis der mechanischen Simulation. Der 
Linsenhalter ist im Betriebsfall (Abbildung 4.6(b)) oben um 10,38 um 
und unten um 6,62 um verschoben. Mit dem Durchmesser des Linsen- 
halters von 10,75 mm ergibt sich eine vernachlässigbare Verkippung 
von 0,35Milligrad. 
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Abbildung 4.5: Links: Isometrische Darstellung der Baugruppe. Rechts: Ansicht von 
vorne. Griine Pfeile: Fixierte Unterseite. 


Eine optische Simulation dieses Ergebnisses wurde nicht durchgeführt, 
da auch die Verschiebung entlang der optischen Achse lediglich ei- 
ne weitere Verringerung der Laserintensität am Konverter nach sich 
zieht. 


Aus den Simulationen folgt, dass in der Montage eines jeden einzel- 
nen TO90-Submounts einmalig die Laserdiode und die Asphäre unter 
Betriebsbedingungen, also unter den im Fahrzeug vorherrschenden 
Bedingungen, bei Betriebsstrom (hier: 2A) und bei Betriebstempe- 
ratur (hier: Tc = 85°C) gegeneinander justiert und fixiert werden 
müssen. 


Diese Erkenntnis führte dazu, dass für die Nachfolgeprojekte der 
im Jahr 2014 in Serie gegangenen u-LARP-Scheinwerfer eine aktive 
Justage von Asphäre gegen Facette in der Fertigung etabliert wurde. 
Dadurch konnten die Leuchtdichte erhöht und die Leuchtdichte- und 
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11.000 iy 
100% 
v 
| 0.0004 + 000 


(a) Gesamtverschiebung bei Te = 25°C (b) Gesamtverschiebung bei Tc = 85°C 
mit 25°C als Referenztemperatur. mit 25 °C als Referenztemperatur. 


Abbildung 4.6: Während die Verschiebung bei Tc = 25 °C noch nur wenige Hundertstel 
pm beträgt, so ist bei Tc = 85°C erkennbar, dass die Oberseite der Linse gegen die 
Unterseite um einige pm verkippt ist. 


Spotpositionstoleranz zweier p-LARP-Gen2-Module verringert wer- 
den. Dadurch lassen sich kleinere Optiken und/oder größere Leucht- 
weiten, erzielen. 


4.3 FOKUSPOSITION UBER VORWÄRTSSTROM 


Für die Justage der Asphäre im Bezug zur Laserdiode wäre es am 
einfachsten, man würde die Laserdiode knapp über dem Schwellstrom 
betreiben. Es ist jedoch wohlbekannt, dass in Halbleiterlasern die 
Halbwertsbreite [11] und die Strahldivergenz [12] und damit auch die 
Fokusposition vom Vorwärtsstrom abhängt. 
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Außerdem ändert sich die Strahldivergenz der fast und der slow axis 
mit der Temperatur [12] und dem Strom [11] unterschiedlich. Auch der 
Thermal Lensing Effekt, der auf eine Erhohung des Brechungsindex 
bei hohen Temperaturen in der aktiven Zone zurtick zu fiihren ist 
[30] hat Auswirkungen auf die Spotposition entlang der optischen 
Achse. 


Daher soll hier messtechnisch untersucht werden, wie groß der Effekt 
für das in diesem Projekt verwendete optische System im statischen 
sowie im dynamischen Fall ist. 


4.3.1 VERSUCHSAUFBAU UND -DURCHFÜHRUNG 


Abbildung 4.7: Messaufbau zur Bestimmung von Leuchtdichte, Temperatur und Leis- 
tungsdichte. 

1: Leuchtdichtekamera / Thermokamera / Ulbrichtkugel; 2: Leistungsmessgerät; 3: Leis- 
tungsmesskopf (Ophir L50(150)A); 4: Ablenkspiegel; 5: Laserdiode mit Asphäre und 
Wasserkühlung; 6: 2-D Scanspiegel; 7: Konverterhalter / CCD-Kamera 
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In Abbildung 4.7 ist die zehnte Iteration des allgemeinen Messaufbaus 
zur Messung der Laserleistung, der Intensitätsverteilung des blauen 
Laserspots, der Leuchtdichteverteilung auf dem Konverter, des Licht- 
flusses nach dem Konverter und der Temperaturverteilung auf dem 
Konverter dargestellt: 


Der Laserstrahl wird von der wassergekühlten Laserdiode, die mit 
einer fokussierenden Asphäre versehen ist, emittiert (5). Dieser wird 
über einen 75/25-Spiegel (4) geteilt. Der 25 %-Anteil trifft auf den Scan- 
spiegel (6), der 75 %-Anteil auf den Leistungsmesskopf (3). Die Mess- 
unsicherheit des Leistungsmesskopfs Ophir L50(150)A beträgt 3% [31]. 
Es folgt der Konverter im Konverterhalter (7). 


An Stelle des Konverterhalters wurde in diesem Fall eine CCD- 
Kamera, BGP-USB-SP620 der Firma Ophir, montiert. Diese kann über 
den x-y-z-Tisch verfahren werden. Die Aufnahme der Kamera ist 
so konstruiert, dass die Chipebene möglichst gut mit der Konver- 
terebene übereinstimmt. Der Verfahrweg in Richtung Laserdiode 
wird mittels einer Messuhr aufgenommen. Position (1) bleibt unbe- 
setzt. 


Auf eine Kompensierung der Junction-Temperatur für unterschiedli- 
che Vorwärtsströme wurde verzichtet, da eine Messung der Junction- 
Temperatur nicht möglich war. Daher sind die Effekte von Vorwärtsstrom- 
und Temperaturänderung überlagert. Die Gehäusetemperatur der 
Laserdiode wurde für alle Messungen gleichbleibend auf 25 °C gehal- 


ten. 


Im statischen Fall wird der Scanspiegel (Position (6) im Abbildung 
4.7) nicht betrieben und der Laserspot ruht in seiner Position auf dem 
CCD-Chip. Für die dynamischen Messungen wird der Scanspiegel 
aktiviert und lässt den zuvor ruhenden Laserspot über den Detektor 
wandern. 
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Zur Auswertung der Laserspotgröße wurde die in der ISO 11146-1: 
2005 als Standard definierte D4r-Methode gewählt. Für einen Gauß- 
förmigen Laserspot entspricht der D40-Wert dem von 1/e? was für 
die Schnitte im dynamischen Fall angewandt wird und die Breite des 
Laserspots an der Stelle angibt, an der die Intensität auf 13,5 % ihres 
Maximalwertes abgefallen ist. 


4.3.2 CHARAKTERISIERUNG DES MESSAUFBAUS 


Um die Leistungsdichte der gemessenen Laserspots berechnen zu kön- 
nen, müssen die von der Kamera in Counts ausgegebenen Werte in 
Watt pro Quadratmeter umgerechnet werden. Dafür muss bekannt 
sein, wieviel Lichtleistung auf den CCD-Chip trifft. Mit dieser Infor- 
mation kann die Software der Kamera die Leistungsdichte berech- 


nen. 


Um den CCD-Sensor der Kamera vor zu hohen Leistungsdichten zu 
schützen, wurde ein Strahlteiler verwendet, der 10 % der Laserinten- 
sität in Richtung der Kamera reflektieren sollte. Die restlichen 90 % 
sollten in einer Strahlfalle deponiert werden. 


Messungen zeigten, dass der Strahlteiler bei der hier verwendeten Wel- 
lenlänge von ca. 445 nm 25 % statt nur 10 % reflektiert. Die gemessenen 
Werte können in Tabelle 4.2 nachgelesen werden. Die Leistungsdif- 
ferenz von Asphäre zu Konverter im 100 %-Fall geht im Scanspiegel 
verloren. Zum Schutz des CCD-Chips mussten weitere Neutraldichte- 
Filter verwendet werden. 
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Tabelle 4.2: Charakterisierung der Leistungsverluste des Aufbaus in Abbildung 4.7. 


Strom Leistung [W] Leistung [W] Leistung [W] Leistung [W] 
[A] nach Asphäre 90/10 Spiegel 100% Spiegel 100% Spiegel 
@ Konverter @ Konv. stat. @ Konv. dyn. 


0,5 0,44 0,09 0,37 0,36 
15 2,03 0,42 1,71 1,70 
2,0 2,73 0,57 2,28 2,28 


4.3.3 STATISCHE MESSUNG 


Zur Berechnung der Leistungsdichte mit dem 100 %-Spiegel wurde 
die gesamte gemessene Leistung aus Tabelle 4.2 auf die in Abbildung 
4.9 gezeigten weißen Rechtecke angewandt. 


Einige der Positionen der Messreihe bei 1,5 A sind in Abbildung 4.9 
zu sehen. Man erkennt oberhalb des Laserspots weitere, sekundäre 
Reflexe aus dem Strahlteiler. 


Abbildung 4.8 zeigt die Breite des Laserspots über den Verfahrweg 
der CCD-Kamera entlang der optischen Achse bei verschiedenen Vor- 
wärtsströmen. Der Fokuspunkt liegt dort, wo die Breite ihr Minimum 
erreicht. Die Position des Minimums variiert mit dem Vorwärtsstrom 
um mehrere Millimeter. Der Verfahrweg des x-y-z-Tisches reichte nicht 
aus, um auch den „fallenden Ast“ bei Ip = 2,0 A aufzunehmen. Die 
Fehlerbalken wurden der besseren Lesbarkeit halber nicht eingetragen. 
Die Ablesegenauigkeit der Messuhr beträgt 5 um, die Standardabwei- 
chung bei der Bestimmung von D4oM ca. 0,1 um, 0,2 um und 0,1 ym 
für 0,5 A, 1,5 A und 2,0 A. Für D4om betragen die Standardabweichun- 
gen entsprechend ca. 0,3 um, 0,4 um und 0,1 ym. 
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Betrachtet man die um mehrere Millimeter unterschiedlichen Fokuspo- 
sitionen für die verschiedenen Vorwärtsströme, lässt sich folgern, dass 
für Anwendungen in denen höchste Leuchtdichten gefordert sind, 
die Asphärenposition genau auf den Vorwärtsstrom angepasst sein 


muss. 


Eine solche Nutzung ist beispielsweise das LARP-Zusatzfernlicht, wie 
es seit 2014 im Markt zu finden ist, das Sichtweiten von über 600 m er- 
möglicht [32]. Weitere LARP-Anwendungen sind Konverterräder, wie 
sie in Projektoren oder zur Endoskopie eingesetzt werden, bei denen 
Leuchtdichten von über 500 cd - mm? oder gar über 1000 cd - mm? 
gefordert sind, was selbst die leistungsstärksten LEDs heute nicht zu 
leisten vermögen [33]. 
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Abbildung 4.8: Breite des Laserspots tiber den Verfahrweg der CCD-Kamera entlang 
der optischen Achse bei verschiedenen Vorwärtsströmen. Der Fokuspunkt liegt dort, 
wo die Breite ihr Minimum erreicht. Die Position des Minimums variiert mit dem 
Vorwärtsstrom. 
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(a) -3,0 mm; D4o,, = 754m; (b) -2,5mm; D4o,, = 60,4um; (c) -2,1mm; D4o,, = 47,9 pm; 
D4oy = 706 pm. D4om = 702 um. D4om = 700 um. 


(d) -1,8mm; D4o,, = 37,8um; (e) -1,5mm; D4¢,, = 31,lum; (f) -1,2mm; D4o,, = 27,1 pm; 
D4om = 696 pm. D4om = 694 pm. D4oy = 693 pm. 


(g) -1,0mm; D4om = 26,0pm; (h) -0,9mm; D4o,, = 25,7 9m; (i) -0,6 mm; D4o,, = 27,0 pm; 
D4om = 694 pm. D4om = 694 ym. D4om = 696 pm. 


(j) -0,2mm; D4o,, = 32,1um; (k) 0,1 mm; D4o,, = 41,lum; (I) 0,4mm; D4o,, = 49,4 am; 
D4oy = 697 um. D4oy = 689 pm. D4oy = 688,6 um. 


(m) 0,7mm; D4o,, = 57,6um; (n) 1,0mm; D4o,, = 67,0 pm; (0) 1,5 mm; D4o,, = 83,0 um; 
D4oy = 688,4 um. D4oy = 687 um. D4oy = 685 um. 


Abbildung 4.9: Statische Laserspotbreitenmessung in unterschiedlichen CCD- 
Positionen entlang der optischen Achse bei Ir = 1500 mA. 
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4.3.4 DYNAMISCHE MESSUNG 


Für die dynamischen Messungen wurde die Scanlänge so eingestellt, 
dass sie möglichst die gesamte Sensorbreite nutzt. Dies ist in Ab- 
bildung 4.10 zu sehen. Zur Umrechnung von Counts in Leistungs- 
dichte wurde die in Tabelle 4.2 angegebene, gemessene optische Leis- 
tung für die verschiedenen Vorwärtsströme gleich der Leistung in- 
nerhalb der in Abbildung 4.10 zu sehenden weißen Rechtecke ge- 
setzt. 


Oberhalb der Scanlinie in Abbildung 4.10 sind, wie bei der statischen 
Messung auch, die zusätzlichen Reflexe des Strahlteilers zu erken- 
nen. 


Aus den Abbildungen 4.10 bis 4.16 lasst sich ermitteln, dass die Leis- 
tungsdichte senkrecht zum Scanweg zwar mit dem Vorwärtsstrom 
variiert, jedoch kaum mit dem Verfahrweg entlang der optischen Ach- 
se. Gleiches gilt für die Scanlinienbreite in 1/e?, gemessen in der Mitte 
der Scanlinie. 


Dies kommt von der leichten Schrägstellung des Laserspots im Be- 
zug zur Scanrichtung. Die slow axis bestimmt die Scanlinienbreite 
und damit auch die Leistungsdichte. Wie in Abbildung 4.8 zu sehen 
ist, variiert die slow axis nur wenig mit dem Verfahrweg entlang der 
optischen Achse. Die Messergebnisse sin in Tabelle 4.3 zusammenge- 
fasst. 
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Tabelle 4.3: Auswertung der Schnitte der dynamischen Leistungsdichtemessung. Die 
Leistungsdichte ist die maximale gemessene Leistungsdiche, die Scanlini- 
enbreite ist der Mittelwert der gemessenen Breiten aller CCD-Positionen 
entlang der optischen Achse, die maximale Abweichung beschreibt die 
Differenz zwischen größter und kleinster gemessener Scanlinienbreite. 


Strom Leistungsdichte Scanlinienbreite 1/e max. Abweichung 


A W- cm~? um um 
0,5 49 169,4 13 
15 201 173,0 18 
2,0 267 171,7 26 
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(a) 0,0 mm. 


(b) -0,5 mm. 


(c) -0,7 mm. 


(d) -0,9 mm. 


(e) -1,1 mm. 


(f) -1,3 mm. 


(g) -1,75 mm. 


Abbildung 4.10: Dynamische Laserspotbreitenmessung mit über den CCD-Chip scan- 
nendem Laserspot in unterschiedlichen CCD-Positionen entlang der optischen Achse 
bei Ip = 1500 mA. 
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Abbildung 4.11: Schnitt durch die Mitte des Scanwegs bei Ir = 0,5 A. Die Leistungs- 
dichte und die Breite senkrecht zur Scanrichtung variieren kaum über den Verfahrweg 


außer in den Extremwerten. 
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Abbildung 4.12: Normalisierter Schnitt durch die Mitte des Scanwegs bei Ir = 0,5 A. 


Die Scanlinienbreite 1/e? variiert wenig über den Verfahrweg. 
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Abbildung 4.13: Schnitt durch die Mitte des Scanwegs bei Ir = 1,5 A. Die Leistungs- 
dichte und die Breite senkrecht zur Scanrichtung variieren kaum über den Verfahrweg. 
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Abbildung 4.14: Normalisierter Schnitt durch die Mitte des Scanwegs bei Ir = 1,5A. 


Die Scanlinienbreite 1/e variiert nur wenig über den Verfahrweg. 
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Abbildung 4.15: Schnitt durch die Mitte des Scanwegs bei Ir = 2,0 A. Die Leistungs- 
dichte und die Breite senkrecht zur Scanrichtung variieren kaum über den Verfahrweg. 
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Abbildung 4.16: Normalisierter Schnitt durch die Mitte des Scanwegs bei Ir = 2,0 A. 


Die Scanlinienbreite 1/2 variiert nur weni über den Verfahrweg. . 
Für Anwendungen in denen höchste Leuchtdichten gefordert sind 
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oder für optische Konzepte mit möglichst kleinen, übereinander ange- 
ordneten Scanlinien muss die Asphärenposition genau auf den Vor- 
wärtsstrom angepasst sein. Außerdem muss dabei darauf geachtet 
werden, die slow axis sehr genau entlang der Scanlinie auszurichten. 
Ist der Laserspot auch nur geringfügig verdreht und wird dann über 
den Konverter gescannt, ist die Scanlinienbreite von der Breite der 


slow axis dominiert. 


Setzt man im Gegensatz dazu die fast axis mit der Scanrichtung gleich, 
ergibt sich ein gegen eine Rotation des Laserspots wesentlich unemp- 
findlicheres System. Da sich die Leistungsdichte innerhalb des Scan- 
wegs kaum über die Verschiebung entlang der optischen Achse ändert, 
könnte ein weiterer Vorteil dieser Anordnung darin liegen, die punk- 
tuelle Leistungsdichte über eine Defokussierung zu reduzieren, um 
optisches Quenching (deutsch: optisches Löschen) des Leuchtstoffs zu 
verringern (siehe Abschnitt 5.3). Dies müsste in weiteren Messungen 


untersucht werden. 


Eine weitere Folge der sich kaum ändernden Leistungsdichte und 
Scanlinienbreite über die Verschiebung entlang der optischen Achse ist, 
dass in der Simulation nur eine geringe Änderung der Temperatur für 
die gemessenen Scanlinien innerhalb einer Vorwärtsstrom-Messreihe 


zu erwarten ist. 


4.4 ZUSAMMENFASSUNG 


Das System aus der hier verwendeten Laserdiode im TO-90 Ge- 
häuse und der Asphäre ist empfindlich gegen kleinste Positions- 
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änderungen der Laserfacette und gegen Änderungen im Vorwärts- 


strom. 


Daher muss in der Montage eines jeden einzelnen TO90-Submounts 
einmalig die Laserdiode und die Asphäre unter Betriebsbedingungen, 
also unter den im Fahrzeug vorherrschenden Bedingungen, bei Be- 
triebsstrom (hier: 2 A) und bei Betriebstemperatur (hier: Tc = 85°C) 


gegeneinander justiert und fixiert werden. 


Die hier gemessenen Leistungsdichten des Laserstrahls werden für die 
Simulation des Konverters als Eingangsgröße im folgenden Kapitel 


Verwendung finden. 
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Kapitel 5 


DYNAMISCHE MODELLIERUNG DES 
LEUCHTSTOFFS 


Die Diffusionsnäherung ist eine unter anderem aus der Biologie be- 
kannte Näherung der Strahlungstransportgleichung. Sie wird zur Be- 
rechnung von Photonentransport in Gewebe genutzt [34]. Auch in 
der Astrophysik findet die Diffusionsnäherung zur Berechnung von 
Strahlungstransport in Sternatmosphären, Nebulae, Novae oder Su- 
pernovae Verwendung [35]. 


Für die Simulation von Strahlungstransport in keramischen Leucht- 
stoffen ist die Diffusionsnäherung erst seit kurzem in der Literatur 
zu finden. Lenef etal. haben sich mit diesem Thema in [13] und [4] 
befasst. In diesen Veröffentlichungen wurde die Diffusionsnäherung 
stets auf einen zweidimensionalen, rotationssymmetrischen Fall ange- 


wandt. 


In dieser Arbeit wird die Simulation keramischer Leuchtstoffe mittels 
der Diffusionsnäherung von rotationssymmetrischen Fällen auf all- 
gemeine, dreidimensionale Problemstellungen erweitert. Außerdem 
wird dieser Leuchtstoff erstmalig mittels der Diffusionsnäherung mit 
einer zeitlich und örtlich veränderlichen Laserlichtquelle simulatorisch 
untersucht. 
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5.1 EINLEITUNG 


Laser Activation of a Remote Phosphor (LARP) ist ein hochgradig 
effektiver Ansatz, sehr hohe Leuchtdichten mit Halbleiterlaserdioden 
zu erzeugen. Diese Lichtquellen haben, verglichen mit LEDs, ein sehr 
viel kleineres Etendue und übertreffen diese in der Leuchtdichte min- 
destens um den Faktor zwei [4]. Dies kann zu einer wesentlich hö- 
heren Leuchtdichte in Projektor- oder Display-Anwendungen füh- 
ren. Auf dem Markt befinden sich heutzutage hauptsächlich LARP- 
Anwendungen, die ein Leuchtstoff-Pulver in einer Silikonmatrix auf 
einem sich drehenden, hoch reflektiven Rad implementieren. Durch 
die Rotation des Rades wird eine sonst punktuelle Strahlenbelastung 
auf den Umfang des Rades verteilt. Dies hilft dabei optische sowie 
thermische Löschung in Grenzen zu halten [4]. 


Um noch höhere Bestrahlungsleistungen erreichen zu können, um sich 
bewegende Teile zu vermeiden und um Bauraum zu sparen, hat OS- 
RAM statische LARP-Systeme, basierend auf einem in hohem Maße 
wärmeleitfähigen, monolithischen, keramischen Leuchtstoff entwi- 
ckelt. Diese Leuchtstoffe sind bereits in einem sehr kompakten und da- 
bei gleichzeitig besonders leistungsfähigen Lichtmodul für Industrie- 
und Produktionsanlagen im Einsatz, das darüber hinaus für mikro- 
skopische und endoskopische Anwendungen genutzt wird. Mit dem 
BMW i8 und dem Audi R8 LMX kamen beinahe zeitgleich die welt- 
weit ersten Serienfahrzeuge mit Laser-Scheinwerfern mit keramischen 
Leuchtstoffen auf den Markt. 


Die optische Gesamtleistung im iLaS Laserscanner übersteigt die der 
sich bereits in Produktion befindlichen Lasermodule um ein Vielfa- 
ches. Daher muss untersucht werden, ob das thermische Verhalten 
des Konvertersystems den sicheren Einsatz im Straßenverkehr ge- 
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währleisten kann. Dazu werden kurz die keramischen Leuchtstoffe 
eingeführt, das sogenannte Quenching erklärt, die Diffusionsnähe- 
rung mit spekularen Randbedingungen eingeführt, deren Einbindung 
in Comsol Multiphysics® aufgezeigt und verschiedene Simulationen 
durchgeführt. 


5.2 KERAMISCHE LEUCHTSTOFFE 


Keramiken sind Materialien, die zumindest aus einem metallischen 
und einem nicht-metallischen Element bestehen. Wenn man es ge- 
nau nimmt, so sind auch die in einer Silikonmatrix eingebrach- 
ten Leuchtstoff-Pulver Keramiken. Hier jedoch soll von kerami- 
schen Leuchtstoffen als polykristalline, dotierte Festkörper die Rede 


sein. 


Spricht man von Ce:YAG oder Ce:GdYAG, so ist zu erwähnen, dass 
diese Materialien kubische Kristallsymmetrie aufweisen und damit 
optisch isotrop sind [4], [36]. Dies impliziert, dass Streuung nicht 
an Korngrenzen stattfindet [4]. Diese findet folglich hauptsächlich 
an Poren statt, die an der Schnittstelle von drei oder mehr Körnern 
entstehen [4]. Sie können demgegenüber auch innerhalb von Körnern 
als Resultat des Sinterprozesses vorkommen [14]. Gelegentlich kann 
es geschehen, dass die in der Keramik verwendeten Chemikalien 
anders als vorgesehen reagieren und innerhalb der Mikrostruktur in 
die gasförmige Phase übergehen [4]. Diese Gaseinschlüsse tragen auch 
zur Streuung bei. Jedoch ist der Brechungsindex im Vergleich zu dem 
der Poren (n = 1) größer und somit die Streuwirkung kleiner, da der 
Unterschied im Brechungsindex geringer ist. 
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Wichtig für diese Arbeit sind die Keramikleuchtstoffe vor allem des- 
halb, weil sie, abhängig von ihrer Zusammensetzung, eine hohe Wär- 
meleitfähigkeit im Bereich von 1 to 100W - m-1.K1 aufweisen [14]. 
Die Wärmeleitfähigkeit von Leuchtstoff-Pulver in Matrixmateriali- 
en hängt im Wesentlichen von dem verwendeten Trägermaterial ab. 
Im Falle von Silikon ist sie im Bereich von 0,1 to 0.3W-m-!-K7! 


[4]. 


Wärmequellen bei der Erzeugung von weißem Licht über Konversion 
von blauem Laserpumplicht sind die nicht strahlenden Übergänge. In 
Abbildung 5.1 sind die vereinfachten Energieniveaus von Ce** darge- 
stellt. Die Verluste pro Photon berechnen sich nach 


BE =he- ( k : : (5.1) 


Akonvertiert Apump 


un, 5d! 


Oa Ce = 


oe 4f! 


Abbildung 5.1: Vereinfachte Darstellung der Energieniveaus von Ce** nach [37] und 
[13]. t und 7,, sind die Zerfallszeiten der strahlenden (r) und nicht strahlenden (nr) 
Übergänge. o,,. die Absorptions- bzw. Emissionsquerschnitte. 
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5.3 QUENCHING MECHANISMEN IN KERAMISCHEN 
LEUCHTSTOFFEN 


Einer der zentralen einschränkenden Faktoren eines Leuchtstoffs in 
einer LARP-Anwendung ist der Verlust der Konversionseffizienz bei 
intensiver Pumpanregung. Die Hauptverluste sind das thermische 
Quenching, was die Quanteneffizienz bei erhöhten Temperaturen ver- 
ringert und das durch die Pumpintensität verursachte optische Quen- 
ching. [4], [13], [14] 


Thermisches Quenching tritt auf, wenn die Zerfallsraten der nicht 
strahlenden Übergänge des oberen Energieniveaus sich der Zer- 
fallsrate des strahlenden Übergangs annähern. Dadurch wird die 
Quanteneffizienz des Leuchtstoffs verringert [13]. Bei Ce-Granaten 
tritt der Beginn des thermischen Quenchings typischerweise im Be- 
reich von 150 °C bis 300 °C auf [38]. Wird Gadolinium hinzugefügt 
oder die Ce-Konzentration erhöht, setzt des thermischen Quen- 
chings bei höheren Temperaturen ein [14]. Bachmanns Daten zeigen, 
dass thermisches Quenching bereits bei 0,33% Ce abnimmt [38]. Es 
kann durch eine modifizierte Arrhenius-Gleichung genähert werden 
[38]: 


t(T) = E (5.2) 
1+n:-C-etsT 
I 
I(T) = — (5.3) 
1+D:.etsT 


Dabei sind t(T) und I(T) die temperaturabhängige Zerfallszeit 
der nicht strahlenden Übergänge und Intensität und Tọ = T(0) 
und Io = I(0). C und D sind Geschwindigkeitskonstanten, E4 die 
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Aktivierungsenergie für diesen Prozess und kg die Boltzmannkonstan- 
te. 


Eine weitere, leicht modifizierte und einfachere Möglichkeit das 
thermische Quenching zu modellieren, ist in [13] und [39] zu fin- 
den: 


= Wp est (5.4) 


Dabei ist T,; die Zerfallsrate der nicht strahlenden Übergänge 
und Wọ eine Geschwindigkeitskonstante. Siehe auch Abbildung 
5:1: 


Das Quenching bei hohen Pumpleistungsdichten verringert ebenfalls 
die Quanteneffizienz. Dabei erhöht allein eine hohe Pumpintensität die 
nicht strahlenden Verluste. Dies sind Mehr-Photonen-Prozesse, die mit 
der Pumpintensität nichtlinear zunehmen. Im Fall von Ce-Granaten 
sind zwei Mechanismen für diese optischen nichtlinearen Verluste 
wahrscheinlich [40], [14]: 


1. Absorption im angeregten Zustand vom oberen Lumineszenz- 
niveau zum Leitungsband und anschließende Übertragung in 
Defektzustände, die nicht strahlend zerstreuen. 


2. Auger-ähnlicher Energietransfer zwischen angeregten benach- 
barten Ce-Ionen. Dabei überträgt ein Ion seine Energie auf den 
Nachbarn. Dadurch geht es in den Grundzustand über und das 
andere wird in das Leitungsband angeregt. Anschließend geht 
dieses in Defektzustände über. 


Folgt man der Literatur, so ist Absorption im angeregten Zustand der 
dominante Prozess für optisches Quenching [14], [40]. 


Für die simulatorischen Betrachtungen des Leuchtstoffs ist es demnach 
essentiell, sowohl das thermische als auch das optische Quenching zu 
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berticksichtigen. Das thermische Quenching wird in der Simulation 
nach Gleichung 5.4 modelliert. Dabei ist nach [13] und [39] EA ~ 
10 000 cm!. Das optische Quenching wird berücksichtigt, indem eine 
angemessene Quanteneffizienz yog für die in Kapitel 4.3 gemessenen 
Leistungsdichten gewählt wird. 


5.4 DIE DIFFUSIONSNÄHERUNG 


Der Photonentransport in keramischen Leuchtstoffen kann nu- 
merisch mit Monte-Carlo-Simulationen oder analytisch mit der 
Strahlungstransportgleichung (STG) modelliert werden. Als Integro- 
Differentialgleichung (s. 5.8) ist sie schwer zu lösen, ohne Näherun- 
gen einzuführen. Eine solche ist die Diffusionsnäherung. Lösungen 
der Diffusionsgleichung für den Photonentransport sind rechne- 
risch effizienter, aber weniger genau als Monte-Carlo-Simulationen. 


[34] 


Bewegen sich Photonen in einem Medium, beschreibt die STG ma- 
thematisch den Energietransfer. Der Strahlungsfluss durch ein Ele- 
ment mit kleiner Fläche im Strahlungsfeld wird durch die Strahl- 
dichte L(7,3,t) charakterisiert. Sie ist definiert als Leistung pro 
Fläche und Raumwinkel. Dabei bezeichnet 7 die Position, § den 
Richtungs-Einheitsvektor und ¢ die Zeit. Einige andere wichtige 
physikalische Größen basieren auf der Definition der Strahldich- 
te: 


Flussdichte: 
O(7,t) = f L(¥,8,t) dQ (5.5) 
Art 
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Stromdichte: 
IH = f SL(7,8,t) dQ (5.6) 
Art 
Bestrahlung: 
+9 
ine / D(F,t) dt (5.7) 


Die STG ist eine Integro-Differentialgleichung, die die Strahldichte 
L(7,8,t) beschreibt. Sie kann mittels Energieerhaltung hergeleitet wer- 
den (siehe z. B. [34]). Die STG besagt, dass ein Lichtstrahl durch Diver- 
genz und Extinktion (bestehend aus Absorption und Streuung vom 
Strahl weg) Energie verliert und Energie von Lichtquellen im Medium 
und Streuung zum Strahl hin gewinnt [34]. Kohärenz, Polarisation 
und Nichtlinearität werden vernachlässigt [34]. Optische Eigenschaf- 
ten wie der Brechungsindex n, der Absorptionskoeffizient fq, der 
Streukoeffizient us und die Streuanisotropie g werden als zeitinvariant 
angesehen, können jedoch räumlich variieren [34]. Die Streuung wird 
als elastisch angenommen [34]. Die STG, auch Boltzmann-Gleichung 
genannt, lautet wie folgt: 


cot = (5.8) 


e Dabei ist c die Lichtgeschwindigkeit im Medium. 
© Ut = Ha + Ms ist der Auslöschungskoeffizient. 


e P(s’,8) die Phasenfunktion, die die Wahrscheinlichkeit darstellt, 
dass Licht mit Ausbreitungsrichtung $ in den Raumwinkel dO 
um 3 gestreut wird. In den meisten Fällen hängt die Phasen- 
funktion nur vom Winkel zwischen dem gestreuten $’ und den 
einfallenden Richtungen ŝ ab, d.h. P(s’,§) = P(8’ - 5) [34]. Die 
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Streuanisotropie kann ausgedriickt werden als 
g = SE SP - 8) dO [34]. 


e S(7,8,t) beschreibt die Lichtquelle. 


In der STG definieren sechs verschiedene unabhängige Variablen die 
Strahlung an jedem räumlichen und zeitlichen Punkt: 7 = (xy,z)", 
Polar- und Azimutalwinkel 0 und ¢ von $ und die Zeit t. Durch ge- 
eignete Annahmen über das Verhalten von Photonen in einem Streu- 
medium kann die Anzahl unabhängiger Variablen reduziert werden. 
Diese Annahmen führen zur Diffusionstheorie und zur Diffusions- 
gleichung für den Photonentransport. Zwei Annahmen erlauben die 
Anwendung der Diffusionstheorie auf die STG [34]: 


1. Es gibt nur sehr wenige Absorptionsereignisse im Vergleich zu 
Streuereignissen. Ebenso treten nach zahlreichen Streuereignis- 
sen wenige Absorptionsereignisse auf und die Strahlung wird 
nahezu isotrop. Diese Annahme wird auch als Richtungsverbrei- 
terung bezeichnet. 


2. In einem hauptsächlich streuenden Medium ist die Zeit für ei- 
ne wesentliche Änderung der Stromdichte viel länger als die 
benötigte Zeit, eine mittlere freie Weglänge eines Transports zu 
durchqueren. Damit ist über eine mittlere freie Weglänge die 
relative Änderung der Stromdichte viel geringer als eins. Diese 
Eigenschaft wird auch als zeitliche Verbreiterung bezeichnet. 


Beide Annahmen erfordern ein Medium mit hoher Albedo, d.h. ein 
Material hoher Streuung [34], was für die in diesem Projekt verwende- 
ten keramischen Leuchtstoffe zutrifft. 


Die Strahldichte kann mit Hilfe der Kugelflächenfunktionen Yn,m ent- 
wickelt werden. In der Diffusionstheorie wird angenommen, dass die 
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Strahlung weitgehend isotrop ist, sodass nur die isotropen und aniso- 
tropen Terme erster Ordnung verwendet werden [34] (dabei sind Lym 
die Entwicklungskoeffizienten): 


L(7,8,t) =), = EN) (5.9) 


Die Strahldichte wird durch vier Terme ausgedrückt: einer für den 
isotropen Term (n=0) und drei für die anisotropen Terme erster Ord- 
nung (n=1). Nutzt man die Eigenschaften der Kugelflachenfunktio- 
nen und die Definitionen der Flussdichte ®(7,t) und der Stromdich- 
te J(7,t), können die Terme nach [34] wie folgt ausgedrückt wer- 
den: 


Loo(?,t)Yo0(8) = = (5.10) 
1 
L Lim(7t)Y1m($) = ai (At) -8 (5.11) 


m=-1 
Damit kann die Strahldichte angenähert werden [34]: 


RSS ON + TEN: (5.12) 


Art 


Fügt man diese Strahldichte in die STG ein, so kann diese in skalarer 
oder in Vektorform dargestellt werden: 


a(t) 


DAN +V - JF) = SAH) (5.13) 
ay (Ft 3 1 a 
ka + (pa + WJC) + 5VOE7L) =0 (5.14) 
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Die Diffusionsnäherung ist eine Näherung für Systeme, in denen die re- 
duzierten Streukoeffizienten viel größer sind als die Absorptionskoeffi- 
zienten [34]. Außerdem muss die Schichtdicke mindestens einige mitt- 
lere freie Weglängen betragen [34]. Auch diese Voraussetzung erfüllen 
die in iLaS verwendeten keramischen Leuchtstoffe. 


Mit der zweiten Annahme der Diffusionsnäherung, dass die rela- 
tive Änderung der Stromdichte über eine mittlere freie Weglänge 
viel geringer als eins ist: 4J/J < 1, reduziert sich die Vektorform 
der Strahlungstransportgleichung (5.14) zum Fick’schen Gesetz 
[13]: 

VO@(7,t) 


5.15 
Ha + (1— g) Ms) eo) 


It) = x 


Eingesetzt in die skalare Form der Strahlungstransportgleichung (5.13) 
bekommt man die Diffusionsgleichung [34]: 


1 d®(7,f) Į 
c ot 


Ug®(7,t) — V [DV®(Zt)] = S(F,t) (5.16) 


Dabei ist D = 3 (Ha + (1-g)us) der Diffusionskoeffizient und 
u, = (1- g)us der reduzierte Streukoeffizient. 


Auffällig ist, dass es in der Diffusionsgleichung keine explizite Abhän- 
gigkeit vom Streukoeffizienten gibt. Stattdessen erscheint im Ausdruck 
für D nur der reduzierte Streukoeffizient. Dies führt zu einer wichtigen 
Beziehung: Bei einer Änderung der Anisotropie des Streumediums 
bleibt die Diffusion gleich, wenn der verringerte Streukoeffizient kon- 
stant ist [34]. 


Angenommen die Intensität M eines Lichtstrahls in Richtung ŝo sei 
Mo(so), so ist nach dem Lambert-Beer’schen Gesetz [13] 
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aM 
— = (noext +a) Mo. (5.17) 
GEN) 


Dabei ist Cext = Osc + Caps der Auslöschungsstreuquerschnitt mit Streu- 
querschnitt dsc und Absorptionsquerschnitt 7,,, der Streuzentren der 
Dichte n. « ist der Absorptionskoeffizient pro Einheitslänge für die 
Ausbreitung im Medium ohne Streuzentren. 


Für eine Geometrie, wie in Abbildung 5.2(a) gezeigt, in der ein Laser- 
strahl als spekulare Lichtquelle von unten auf das Medium trifft, wird 
es Rückreflexionen geben. Daher müssen beide spekularen Intensitä- 


ten M+ unter Verwendung geeigneter spekularer Randbedingungen 
gelöst werden. 


Mit [41] lässt sich nach [13] zeigen, dass sich die Diffusionsgleichung 
für einen eingehenden, spekularen Lichtstrahl wie folgt modifizieren 
lässt: 


-V(DV®)+w® =e+e5 +. -V (Jp +h)- (5.18) 


Dabei ist u. = Naps + &. € ist die Leistungsdichte einer Volumenquelle. 
Für die Lumineszenz (das konvertierte, emittierte Licht) ist ¢¢(we) 
bei der Frequenz w., angeregt durch die Pumplichtquelle mit der 
Frequenz wp gegeben durch [13] 


w 
Ee(We) = a (Wp )YQE—— (®+M++M_). (5.19) 

p 
Der Anteil des spekularen Lichts, der gestreut wird und diffuses Licht 


erzeugt, wird mittels einer effektiven Quellfunktion eo beschrieben 
[13]: 


£g = Maıno;.. (5.20) 
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Die Stromdichte in der Drift-Diffusionsform ist [13]: 
J =-DV®+ Jo + Jo (5.21) 


mit 


JÈ = +83gDe}. (5.22) 


Vorausgesetzt, dass an den Grenzflächen Energieerhaltung gelten 
muss, kann nach [13] und [42] die Bestrahlungsstärke an der Grenzflä- 
che gleich dem Integral der reflektierten Strahldichte gesetzt werden. 
Für eine Oberfläche mit Reflektionskoeffizient R(n - 5) und Normalen- 
vektor fi gilt die Randbedingung: 


lol Re) 4 taga ER) =0. (5.23) 
mit [42] 
m/2 
Ro = T 2 sin 0 cos ØR presne (0) d0 (5.24) 
0 
n/2 
Ry = f 3 sin 0 cos? OR Fresnel (0) dO (5.25) 


mit dem Fresnel-Reflektionskoeffizient Rpyesne| (0) und 
cos = §- 7. 


Die im Konverter (Index c) und im Kleber (Index g) erzeugte ther- 
mische Verlustleistung kann nun mittels der temperaturabhängigen 
Quanteneffizienz 70g(T), der Pump- und Emissions-Intensitat P, ; ¢ 


und M= ausgedrückt werden: 
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Qth,c = (1 = nor(T)) *Ha,c,p° (®,,c +M? +Mz) + Ha,c,e* Pe,c 


(5.26) 

Qrh,¢ = Ha,g,p° (On + Mg A: Mz) + Ha,g,e Pe, g . (5.27) 
Dabei ist 

nor(T) poe (5.28) 


= 1+7- VW tar (T) 
mit 1/7,,(T) aus Gleichung (5.4). 


Die STG als Integro-Differentialgleichung wurde durch die Diffusions- 
näherung in eine partielle Differentialgleichung überführt. Nun kann 
die Diffusionsgleichung mit spekularen Randbedingungen mittels ei- 


ner Finiten-Elemente-Simulation simuliert werden. 


5.5 SIMULATION IN COMSOL MULTIPHYSICS® 


Comsol Multiphysics® ist ein Simulationsprogramm, das es ermög- 
licht verschiedene physikalische Probleme gleichzeitig zu simulieren. 
In diesem Fall den Wärmetransport und die Lösung der Diffusionsnä- 
herung der Strahlungstransportgleichung. 


5.5.1 AUFBAU DES MODELLS UND RANDBEDINGUNGEN 


Die Simulationen wurden mit Comsol Multiphysics® in der Version 
5.3 durchgeführt. Das Modell ist in Abbildung 5.2(a) dargestellt und 
besteht aus drei Festkörpern: 
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1. dem Trägersubstrat Saphir (Dicke 500 um), 
2. einem Klebstoff (Dicke 15 um, Messung in Abschnitt 5.5.2) 
3. und dem Konverter (Dicke 70 pm). 


Da der Konverter mittels Warmeleitpaste in einen Halter aus Alumi- 
nium eingeklebt ist, sind die Flächen, an denen der Saphir mit dem 
Halter über den Wärmeleitkleber in Kontakt steht mit einem Wär- 
meübergangskoeffizienten von 10000W - m-1.K”? versehen (siehe 
Abbildung 5.2(b) und Abschnitt 2.1.1). Die Fläche hat eine Temperatur 
von Tas = 70°C, da dies die in diesem Projekt erwartete Temperatur 
im Scheinwerfer ist. 


VEN 


(a) Der Simulationsaufbau von oben nach (b) Kontaktfläche des Wärmeleitklebers 
unten: Konverter, Kleber und Saphir. mit dem Saphir. 


Abbildung 5.2: Simulationsaufbau in Comsol Multiphysics®. 


Die Diffusionsgleichung soll in Comsol Multiphysics® für zwei Licht- 
farben mit unterschiedlichen optischen Eingenschaften in den Medien 
Klebstoff und Konverter gelöst werden: 


1. die blaue Laser-Pump-Quelle mit Index p und 
2. das emittierte, konvertierte Licht mit Index e. 


Daher muss je eine „Physik“ dafür angelegt werden. Entsprechend 
der Form der Diffusionsgleichung wählt man die „Coefficient Form 
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PDE”. Dabei steht PDE für partial differential equation oder auf deutsch: 
partielle Differentialgleichung. Comsol Multiphysics® nimmt die PDE 
in der Form von Gleichung 5.29 auf. 


I, je | dye 
ea apr rag F t v-( CV @ y/o — A® y/o +7) + B-V®p5/e 
+ ape = f 
(5.29) 


Mit dem vorherigen Kapitel ergibt sich, dass 
e e, = 0 und a = 0, wobei 0 der Nullvektor ist, 


e ye =0, Yp = (0,0,Jo)', wobei 2 der Einfallsrichtung des Laser- 
strahls entspricht (siehe Abbildung 5.2(a)), 


© d, = 1/cmitc der Lichtgeschwindigkeit, 

e c der Diffusionskoeffizient Dp/e = 3 (pa p/e + (1 = 8)Hs sie); 
è a= Wa, p/e und 

e f = £p/e, wobei im Kleber offensichtlich ee = 0 gilt. 


Die Parameter sind material- und wellenlängenabhängig, sodass sie 
in Comsol Multi-physics® für jede Domäne und jede „Physik“ separat 
angelegt und berechnet werden müssen. Sie sind in der zentralen For- 
schungsabteilung von OSRAM gemessen und an den entsprechenden 
Stellen im Modell eingesetzt worden. 


Die Diffusionsgleichung wird im Kleber und im Konverter gelöst, da 
der Saphir, bis auf die Fresnel-Reflexion am Eintritt, bei dem hier an- 
genommenen, senkrechten Lichteinfall keinen Einfluss auf den Strahl 
nimmt. Außerdem ist er transparent und erfüllt damit nicht die Rand- 
bedingungen für die Diffusionsnäherung. 
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Die Randbedingung aus Gleichung (5.23) ist mit ihren entsprechenden 
Brechungsindizes auf die obere Flache des Konverters und die unte- 
re Flache des Klebers angewandt. Der Einfachheit halber haben die 
dünnen Randflächen des Klebers und des Konverters eine „null Fluss” 
Randbedingung: 


ñ. (-cvo, Je — 4D yje + 7) =0. (5.30) 


Die Anfangswerte für ®, Je und 9®,,./at sind Null. 


Die Wärmeleitung wird mittels der „Physik“ Wéarmeleitung im Fest- 
körper in Comsol Multiphysics® umgesetzt. Die zu lösende PDE ist 
gegeben durch 


CP + eCpil- VT+V-(-AVT)=Q. (5.31) 


Wobei A die Wärmeleitfähigkeit, ọ die Dichte und Cp die spezifische 
Wärmekapazität des jeweiligen Materials darstellen. Q ist für den 
Konverter gegeben durch Gleichung (5.26) und für den Kleber durch 
Gleichung (5.27). 


Als Anfangswert wurde für die Temperatur aller Körper ein Wert 
von To = Tys = 70°C angenommen, da dies die in diesem Projekt 
erwartete Temperatur im Scheinwerfer ist. 


Für die Auslegung des Netzes für den Kleber und den Konverter wur- 
de die in Abschnitt 4.3 aufgezeigte Laserspotgröße für die x- und 
y-Richtung berücksichtigt. Da der Saphir ein transparentes Medi- 
um ist und damit die Annahmen, die zur Diffusionsnäherung ge- 
führt haben nicht angewandt werden können und die Effekte auf 
die hindurchtretende Strahlung vernachlässigbar sind, wurde dort 
nur die Wärmeleitung gelöst. Daher wurden dort größere Zellen ver- 
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wendet, um Rechenzeit zu sparen. Das Netz ist in Abbildung 5.3 
gezeigt. 


(a) Vernetze Baugruppe. (b) Zellgröße im roten Bereich. 


Abbildung 5.3: Das Netz spiegelt mit seinen an den Grenzflächen kleiner werdenden 
Zellen die Physik der Diffusionsnäherung wider. Mit Zellgröße I < 25 ym in x-y. 


5.5.2 SCHICHTDICKENMESSUNG KLEBER 


Da der Klebeprozess von Konverter auf Saphir manuell durchgeführt 
wurde, wurden zur Ermittlung der Kleberschichtdicke Schnittbilder 
mit dem Rasterelektronenmikroskop in der zentralen Forschungsabtei- 
lung bei OSRAM angefertigt. Im Schnitt ergibt sich eine Kleberdicke 
von tg = (15 + 5) um, wobei die Extremwerte zwischen 3 ym und 


30 pm liegen. Diese sind in Abbildung 5.4 zu sehen. 


5.5.3 ZEITLICH AUFGELOSTE SIMULATIONEN 


Zur erstmaligen Simulation eines nicht rotationssymmetrischen, drei- 
dimensionalen, keramischen Leuchtstoffs mit über der Zeit örtlich ver- 
änderlicher Pumpanregung wurden die in Abbildung 5.6 gezeigten La- 
serspots verwendet. Dies war zum Zeitpunkt der Erstellung der Simu- 
lation die angedachte Leistungsdichteverteilung. Alle sechs Laserspots 
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azum goun 30,9um 25,7um 


Abbildung 5.4: Extremwerte der Schichtdickenmessung des Klebers: Links 3 pm und 
rechts 30 um. Man beachte die unterschiedlichen Skalen. 


tragen dieselbe Leistung von 1,4 W, wobei drei eine Flache von As = 
450 - 455 um? bedecken und damit eine größere Leistungsdichte auf- 
weisen als die größeren Laserspots mit A; = 450 - 1006 um?. Scanrich- 
tung ist die in Abbildung 5.6 erkennbare x-Achse. 


Diese Scanrichtung entlang der slow axis stellt besondere Heraus- 
forderungen an die Genauigkeit der Rotation der einzelnen Spots 
gegeneinander (siehe dazu Abschnitt 4.3) und ist damit aus System- 
sicht betrachtet nicht optimal gewählt. Für die Simulation ist dies aber 


irrelevant. 


Die Laserspots werden mit einer Frequenz von 133 Hz mittels einer 
Dreiecksfunktion über den Konverter verfahren (siehe Abbildung 5.5). 
Auch die spätere Applikation arbeitet bei 133 Hz. 
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I 
0 10 20 30 40 50 60 70 80 
Zeit in ms 


Abbildung 5.5: Dreiecksfunktion als Scanfunktion fiir die Laserspots. 0: Laser in der 


Mitte des Konverters, 1/-1: Laser im Umkehrpunkt 1/2. 
2 W/m2 
| x10? m x10° 
6 
5 
| 
4 
3 
2 
1 
5 0 


Abbildung 5.6: Blauer Input für die zeitabhängige Simulation. Sechs Laserspots als 
Tophat mit gleicher Leistung aber unterschiedlicher Fläche: As = 450 : 455 um?, Ay = 
450 - 1006 um?. Scanrichtung ist entlang der x-Achse. 
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Es wurden zwei Simulationen durchgeführt. In der ersten Simulation, 
die die Zeiten 0-9 ms beinhaltet, wurde es Comsol Multiphysics® 
tiberlassen den Zeitschritt automatisch zu bestimmen. Er betrug 
2,74 us. Die Simulation der ersten 9ms dauerte 1557 004s. Das ent- 
spricht rund 18 Tagen und 30 Minuten. Es wurden 3563 Zeitschritte 
berechnet, 190.07 GB an Arbeitsspeicher und 394.67 GB an virtuellem 
Speicher verbraucht. 


Es ist wichtig am Anfang einer zeitabhängigen Simulation kleine Zeit- 
schritte zu machen, um den Fehler, der sich aus der Simulation ergibt 
und sich in spätere Zeiten fortpflanzt, so gering wie möglich zu halten. 
Um einen sanften Anlauf und damit die Stabilität der Simulation zu ga- 
rantieren, wurde der in Gleichung 5.29 zu findende Dämpfungsfaktor 
d, auf eins gesetzt. 


Nachdem die erste Simulation nach 18 Tagen abgebrochen wurde, 
ist in der zweiten Berechnung der Zeitschritt manuell auf 37,59 us 
festgesetzt worden. Dies entspricht 200 Berechnungsschritten pro Peri- 
odendauer, bzw. 92 um überfahrener Konverterstrecke pro Zeitschritt. 
Dies entspricht rund einem Fünftel der Laserspotbreite und gibt damit 
eine ausreichend gute Auflösung nach der anfänglich wesentlich fei- 
neren Granularität der Zeitschritte. Die Berechnung wurde bei 7,5 ms 
fortgesetzt und lief bis 75ms und damit ca. zehn Periodendauern. Der 
leichte Überlapp der beiden Berechnungen wurde genutzt, um zu 
überprüfen, ob sich Abweichungen in den Simulationsergebnissen zei- 
gen, was nach Abbildungen 5.12 bis 5.15 nicht der Fall ist. Simulation 
2 lief 804934 s, was 9 Tagen, 7 Stunden und 35 Minuten entspricht. Sie 
brauchte mit 252.,53 GB fast den gesamten Arbeitsspeicher von 256 GB 
auf und belegte 397,11 GB virtuellen Speichers. 


Die Ergebnisse der Simulationen sind in Abbildungen 5.7 bis 5.16 ge- 
zeigt. Grafiken 5.7 bis 5.11 zeigen die Temperaturverteilungen einiger 
Zeitschritte. Das erste Bild der Serie zeigt das Temperaturfeld auf der 
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Oberfläche des Konverters im ersten Umkehrpunkt. Man erkennt, dass 
die Temperatur zum Rand hin am höchsten ist, was damit zu erklären 
ist, dass die Wärme, die zu Beginn des Scans in der Mitte des Kon- 
verters eingebracht wurde, bereits ein wenig weiter verteilt ist. Das 
zweite Bild der Serie, Abbildung 5.8, zeigt die Temperatur der Konver- 
teroberfläche nach knapp einer halben Periodendauer. Man erkennt 
wie zuvor, dass dort, wo der Laserspot den Konverter gerade überstri- 
chen hat, die Temperatur am höchsten ist. Im dritten Bild der Serie ist 
der Zustand im zweiten Umkehrpunkt zu sehen. Es ist zu erkennen, 
dass die maximale Temperatur im Vergleich zur vorherigen Abbildung 
abgenommen hat, da die Wärme mittels Konvektion abtransportiert 
wurde. Das vierte Bild zeigt die Temperaturverteilung nach knapp 
einer Periode. Man erkennt, dass die Temperatur fast wieder den Wert 
aus Abbildung 5.8 erreicht hat. Das letzte Bild der Serie zeigt die Tem- 
peraturverteilung am Ende der Simulation bei 75 ms. Man erkennt zu 
den Randern der Temperaturverteilung ein auswaschen, was mit dem 
Abtransport der im Konverter und im Kleber erzeugten Wärme zu 
erklären ist. 


Abbildung 5.7: Temperaturverteilung bei 1,875 ms. Erster Umkehrpunkt. 
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Abbildung 5.8: Temperaturverteilung bei 3,75 ms. Eine halbe Periode. 


Abbildung 5.9: Temperaturverteilung bei 5,625 ms. Zweiter Umkehrpunkt. 


Abbildung 5.10: Temperaturverteilung bei 7,5 ms. Eine volle Periode. 
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Abbildung 5.11: Temperaturverteilung bei 75 ms. Zehn Perioden. In y-Richtung erkennt 
man einen ersten Abfall der Temperaturen zum Rand hin. 


Abbildung 5.12 zeigt die maximale Konvertertemperatur, egal wo im 
Konverter diese auftritt. Der Graph korreliert mit den im vorigen Ab- 
satz getätigten Aussagen. Die maximale Temperatur steigt bis zum 
ersten Umkehrpunkt stetig an und nimmt dann sprunghaft zu, da 
die Zeit nicht ausgereicht hat, die entstandene Wärme von der Be- 
wegung in x-Richtung ab zu transportieren und damit der bereits 
erwärmte Konverter weiter aufgeheizt wird. Die maximale Konverter- 
temperatur nimmt darauf folgend wieder ab, da die Zeit des letzten 
Überstreichens immer weiter zurück liegt und damit schon mehr Wär- 
me abtransportiert werden konnte. Der Knick bei 3,76 ms, einer halben 
Periodendauer, zeigt an, dass der Konverter hier noch nicht vom Laser 
überstrichen wurde. Man erkennt, dass die maximale Konvertertem- 
peratur lokale Maxima in den Umkehrpunkten und lokale Minima 
im Nulldurchgang hat. Außerdem ist erkennbar, dass nach 75 ms der 
eingeschwungene Zustand noch nicht erreicht ist. 
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Abbildung 5.13: Simulationsergebnis Verlustleistung Qrp, c+g über t. 
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In Abbildung 5.13 ist die im Kleber und im Konverter erzeugte Ver- 
lustleistung über der Zeit aufgetragen (siehe Gleichungen (5.26) und 
(5.27)). Man erkennt den sanften Anstieg zu Beginn, der auf den zur 
Stabilisierung eingeführten Dämpfungsfaktor d, in Gleichung (5.29) 
zurückzuführen ist. Diesem Dämpfungsfaktor sind auch die Einbrü- 
che in den Umkehrpunkten geschuldet, die nicht der tatsächlichen 
Physik entsprechen. 


Gleiches gilt für Abbildung 5.14, in der die absorbierte Pumpleistung 
über der Zeit dargestellt ist. Es ist zu erkennen, dass die absorbierte 
Pumpleistung und die Verlustleistung korrelieren, was bedeutet, dass 
die Absorption des emittierten Lichts nur einen geringen Beitrag zur 
gesamten Verlustleistung leistet. 


ol 
T 


Absorb. Pumpleistung in W 
w Ss 
T T 


Simulation 1 


Simulation 1 fortgesetzt 


i tatahi Eidal fet atl fe isl 
1071 10° 101 102 
Zeit in ms 


Abbildung 5.14: Simulationsergebnis absorbierte Pumpleistung tiber t. 
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Der Lichtfluss kann näherungsweise bestimmt werden, wenn man 
von einem Lambert’schen Strahler ausgeht. Dazu muss man das Ober- 
flachenintegral über die Konverteroberfläche über die austretende 
Strahlung bilden, welche zuvor mit der Hellempfindlichkeitskurve 
V(A) gefaltet wurde, welche die spektrale Helligkeits-Empfindlichkeit 
des menschlichen Auges bei Tageslicht beschreibt. 


Der normalisierte Lichtfluss in Abbildung 5.15 zeigt das gleiche unphy- 
sikalische Verhalten wie die Verlustleistung und die absorbierte Pump- 
leistung. Es ist erkennbar, dass der normalisierte Lichtfluss, wie erwar- 
tet, mit der absorbierten Pumpleistung korreliert. 


Normalisierter Lichtfluss 


Simulation 1 fortgesetzt | _| 


ot Simulation 1 


1071 10° 101 102 
Zeit in ms 


Abbildung 5.15: Simulationsergebnis normalisierter Lichtfluss tiber t. 
Abbildung 5.16 zeigt alle zuvor diskutieren Größen, normalisiert auf 1. 


Zur Orientierung sind Hilfslinien in den Umkehrpunkten und eine Re- 
präsentation der Laserspotposition eingezeichnet. Man erkennt, dass 
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die Verlustleistung und die absorbierte Pumpleistung aufeinander 
liegen und die Absorption des emittierten Lichts nur einen geringen 
Beitrag zur gesamten Verlustleistung leistet. Der Lichtfluss korreliert 
mit der absorbierten Pumpleistung. Der langsame Anstieg des Licht- 
flusses und die Einbrüche der absorbierten Pumpleistung sind mit 
der zur Stabilisierung eingeführten Dämpfung d, zu erklären. Der 
Lichtfluss und die Verlustleistung verändern sich nicht temperatur- 
abhängig, da die Temperatur weit unterhalb der Aktivierungsenergie 
aus Abschnitt 5.3 liegt. 


—— Temperatur —— Absorbierte Pumpleistung - - - Laserposition 
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roy de int TTTTTTY Tree 
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H 
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= 
H 
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Abbildung 5.16: Simulationsergebnis normalisierter Lichtfluss über t. Laserposition: 
y =0,5: Mitte des Konverters, y=1: Umkehrpunkt 1, y=0: Umkehrpunkt 2. Man erkennt 
die Temperaturabhängigkeit des Lichtflusses und der Verlustleistung. 
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Es wurden weitere Simulationen durchgeführt, in denen die Damp- 
fung d, ausgeschaltet wurde. Die Berechnung wurde mit einem gröbe- 
ren Zeitschritt von 75,12 us durchgeführt, was 100 Berechnungsschrit- 
ten pro Periodendauer, bzw. 184jım überfahrener Konverterstrecke 
entspricht. Die Simulation ergibt teilweise negative Werte für die Pum- 
pintensität &,, was die Dämpfung zuvor verhindert hatte. Dabei ist es 
unerheblich, ob die Simulation bei t=0s gestartet oder bei t=7,5 ms 
fortgesetzt wird. Es muss weiter eruiert werden, wie eine Simulation 
ohne Dämpfung durchgeführt werden kann. 


Mittels des thermischen Widerstands, der aus Abbildung 5.27(a) zu 
1,7 K : W7! berechnet werden kann, der spezifischen Wärmekapazität 
des Saphirsubstrats und dessen Masse kann nach Gleichung (5.32) 
die Zeit abgeschätzt werden, die zum Aufwärmen benötigt wird: 
4t = 33s. Dabei entspricht 4T einer Erwärmung auf 98 %. 3,3 ent- 
sprechen rund 440 Periodendauern. Die Simulation von zehn Perioden- 
dauern hat vier Wochen benötigt. Mit zukünftigen Rechenleistungen 
könnte es möglich sein, diese langen Berechnungsdauern zu verkür- 


zen. 
T=Ry-Cy-m (5.32) 


ZUSAMMENFASSUNG 


Es wurde zum ersten Mal die Anwendung der Diffusionsnäherung 
der Strahlungstransportgleichung auf einen dreidimensionalen, nicht 
rotationssymmetrischen keramischen Leuchtstoff mit über der Zeit 
örtlich veränderlicher Laseranregung gezeigt. Die analysierten Größen 
verhalten sich wie erwartet. Es bleibt zu klären, wie eine zeitabhängige 
Simulation ohne Dämpfung realisiert werden kann. 
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Zusätzlich zu einem scannenden Lasersystem lässt sich die zeitab- 
hängige Simulation mittels der Diffusionsnäherung auch auf Phos- 
phorräder, wie sie in Projektoren vorkommen oder in der Endoskopie 
Verwendung finden, anwenden. Dazu müssen die entsprechenden 
Streukoeffizienten und Reflexivitäten gemessen werden, sollte es sich 
um in Silikon oder Wasserglas eingebettete Konverter handeln. Diese 
müssen darauf folgend auf ihre Tauglichkeit für die Diffusionsnä- 
herung untersucht und entsprechend ins Modell eingesetzt werden. 
Hier lassen sich eventuell Symmetrieeigenschaften nutzen, um die 
Berechnungszeit zu verkürzen. 


Ein weiteres Anwendungsgebiet könnte die Analyse von keramischen 
Konvertern auf LEDs im pulsweitenmodulierten Betrieb sein. Da hier 
sehr ähnliche Konverter Verwendung finden, ist lediglich die Rand- 
bedingung der spekularen Pumpanregung zu ändern. In [13] wird 
dagegen argumentiert, dass die spekularen Randbedingungen beibe- 
halten werden können, da der spekulare Lichteinfall in einem hoch 
streuenden Medium, wie dem keramischen Konverter, schnell diffus 
wird, sodass der Einfluss der Randbedingung als nicht signifikant 
betrachtet werden kann. Auch bei der Simulation von LEDs lassen 
sich die Symmetrieeigenschaften der LED nutzen, um Simulationszeit 


zu sparen. 


Betrachtet man auch Alterungsmechanismen des Konverters und passt 
die Streukoeffizienten in der Simulation entsprechend an, kann man 
auch gealterte, hoch streuende, keramische Konverter, egal in welcher 
Anwendung, simulieren, und die Unterschiede zwischen Tag eins und 
am Ende der Lebensdauer herausarbeiten. 
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5.5.4 ZEITLICH GEMITTELTE SIMULATIONEN 


Da die Berechnungsdauer der zeitabhängigen Simulationen für zehn 
Periodendauern vier Wochen betragen hat, wurde eine schnellere 
Alternative gesucht. Diese wurde in der über die Zeit gemittelten dy- 
namischen Leistungsdichte gefunden. Auch dies stellt eine Neuerung 
dar, da bisher noch nie dreidimensionale, nicht rotationssymmetri- 
sche, keramische Konverter mittels der Diffusionsnäherung simuliert 
wurden. Die Simulationen dauern im Schnitt 25 Minuten und brau- 
chen mit rund 140 GB Arbeitsspeicher und 160 GB virtuellem Speicher 
wesentlich weniger Ressourcen. 


Abbildung 5.17 zeigt eine der in Abschnitt 4.3 gezeigten, gemessenen 
Scanlinien. Für die Simulation wurden einige dieser Leistungsdichte- 


verteilungen in Comsol Multiphysics® 


eingelesen. Da in den Ergebnis- 
sen der Simulationen für 1,5 A, die zuerst durchgeführt wurden, die 
Temperatur nur um 0,1°C variierte, sind für die Simulationen bei 0,5 A 
und 2,0 A nur wenige Stützstellen simuliert worden. Dieses Verhalten 
war nach der Diskussion der Messergebnisse in Abschnitt 4.3.4 nicht 
anders erwartet worden, musste aber dennoch überprüft werden. In 
den Simulationen bei 0,5 A und 2,0A ergaben sich Abweichungen 
von unter 0,2°C. Die Rückseitentemperatur des Konverters wurde 
auf 65°C geändert, da sich die Anforderungen im Projekt geändert 


hatten. 


Die Grafiken in Abbildung 5.18 zeigen die Temperaturverteilungen 
auf der Konverteroberfläche bei den verschiedenen Vorwärtsströmen, 
wobei die Temperaturskala das jeweilige Maximum im Plot darstellt. 
Für diese Darstellungen wurde für den Import der Laserleistungsdich- 
te wurde die Position entlang der optischen Achse gewählt, die in 
Abschnitt 4.3 als Fokusposition ermittelt worden war. Man erkennt, 
dass sich die relative Temperaturverteilung nicht ändert. Auch das 
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Abbildung 5.17: Laserinput in W - m~?, 2000 mA, 2280 mW, Pos: —2900 pm. Beispiel- 
hafter dreidimensionaler Plot einer Messung aus Abschnitt 4.3. 


liegt an der in Abschnitt 5.17 gezeigten, sich kaum verändernden Sc- 
anlinienbreite. Aus demselben Grund sind keine weiteren Ergebnisse 
weiterer Stützstellen aufgezeigt. 


In Abbildung 5.19 sind dieselben Bilder wie zuvor auf die maximale 
Temperatur in 5.19c skaliert. Es ist festzustellen, dass die Temperatur, 
wie erwartet, mit dem Vorwärtsstrom steigt. 


Damit kann man berechnen, dass der thermische Widerstand kaum 
mit der eingebrachten Leistung variiert (siehe Tabelle 5.1). Dies ist 
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erwartet worden, da die thermischen Leitfähigkeiten innerhalb des 
relevanten Temperaturbereichs als konstant angenommen wurden 
und sich die Fläche unter der Scanlinie nicht signifikant ändert. Siehe 
dazu auch Gleichung (2.1) zur Berechnung des thermischen Wider- 
standes. 


Tabelle 5.1: Ergebnisse der Simulation mit zeitlich gemitteltem blauen Input. Die ma- 
ximale Konvertertemperatur variiert mit dem Vorwärtsstrom, wobei der 
thermische Widerstand weitestgehend konstant bleibt. 


IF [A] PLaser [W] Tmax [°C] Rih [K : wl] 


0,5 0,37 68,1 84 
15 1,70 79,1 8,3 
2,0 2,28 83,8 8,3 
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(a) 0,5 A, 370 mW, Pos: —1300 pm. 


(b) 1,5 A, 1700 mW, Pos: —900 pm. 


(c) 2,0 A, 2280 mW, Pos: —2900 pm. 


Abbildung 5.18: Ergebnisse der zeitlich gemittelten Konvertersimulation. Temperatur- 
skala jeweils auf Maximum. Die relative Temperaturverteilung ändert sich nicht. 
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(a) 0,5 A, 370 mW, Pos: —1300 um, Tmax = 68,1 °C. 


183.3 
rC 


(c) 2,0 A, 2280 mW, Pos: —2900 um, Tmax = 83,8 °C. 


Abbildung 5.19: Ergebnisse der zeitlich gemittelten Konvertersimulation. Temperatur- 
skala jeweils auf 85°C Maximum. Als Position wurde die in Abschnitt 4.3 ermittelte 
Fokusposition gewählt. Die maximale Temperatur ändert sich mit der eingebrachten 
Laserleistung. 
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THERMISCHES QUENCHING 


Um das Phänomen des thermischen Quenchings zu zeigen, wurden 
Simulationen bei gleichbleibenden Randbedingungen aber verschiede- 
nen Kühlflächentemperaturen durchgeführt. 


Auch für die näherungsweise Berechnung der Farbe kann von einem 
Lambert’schen Strahler ausgegangen werden. Damit kann man mittels 
des Oberflächenintegrals der aus der Konverteroberseite austretenden 
Strahlung nach DIN EN ISO/CIE 11664-3:2020-03 [43] die X, Y, Z Werte 
der emittierten Strahlung ermitteln. Der Einfluss der Temperatur auf 
die Farbe ist in Abbildung 5.20 verdeutlicht: Wird der Konverter heiß, 
verringert sich die Quanteneffizienz der Konversion und damit der 
Anteil an konvertiertem Licht, was neben eins Abfalls des Lichtflusses 
auch eine Farbveränderung hervorruft. Dies kann soweit gehen, dass 
die gesetzlichen Vorschriften für die Lichtfarbe „weiß“ (siehe [44]) 
nicht mehr erfüllt sind. 
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Abbildung 5.20: Farbe x und y über der maximalen Konvertertemperatur Tax bei 
Ip = 1,5A. Man erkennt das thermische Quenching und den damit verbundenen 
Farbdrift bei höheren Temperaturen. 
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03 04 05 0.6 


Abbildung 5.21: Temperaturserie der Farbe x und y bei Ir = 1,5 A im CIE Farbdia- 
gramm. Gezeigt sind die Farbmesspunkte und die maximale Anfangs- und Endtempera- 
tur Tmax des Konverters aus der Messserie. Des weiteren ist das in der ECE Regelung Nr. 
48 [44] in Sektion 2.29.1. definierte Weißfeld eingezeichnet. Man erkennt das thermische 
Quenching und den damit verbundenen Farbdrift ins Blaue bei höheren Temperaturen. 
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Bei konstanter Pumpleistung und sich erhöhender die Konverter- 
Kühlflächen-Temperatur Txr, nimmt die Differenz aus der maximalen 
Konvertertemperatur und Tgr in Abbildung 5.22 zu. 


Dasselbe Verhalten zeigt die im Konverter erzeugte Verlustleistung 
in Abbildung 5.23. Der normalisierte Lichtstrom in Abbildung 5.24 
nimmt mit steigender Tgr ab. Dies ist mit der nach Gleichung (5.28) 
verringerten Quanteneffizienz bei hohen Temperaturen zu erklären. 
Man erkennt, dass nach Möglichkeit eine maximale Konvertertem- 
peratur von 150 °C nicht überschritten werden sollte, wenn man den 
Konverter nach dem Modell aus Gleichungen (5.4) und (5.28) nicht zu 
ineffizient betreiben will. 
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Abbildung 5.22: Temperaturdifferenz der maximalen Konvertertemperatur und der 
Kühlflächentemperatur Tmax — Txp über Kühlflächentemperatur Tgr bei Ip = 1,5A. 
Man erkennt das thermische Quenching bei höheren Temperaturen. 
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Abbildung 5.23: Verlustleistung über der maximalen Konvertertemperatur Tmax bei 
Ip = 1,5A. Man erkennt das thermische Quenching bei höheren Temperaturen. 
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Abbildung 5.24: Normalisierter Lichtstrom über der maximalen Konvertertemperatur 
Tmax bei Ir = 1,5 A. Man erkennt das thermische Quenching bei höheren Temperaturen. 
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VARIATION DER KLEBERDICKE 


In Abschnitt 5.5.2 wurde ermittelt, dass die Kleberdicke tg zwischen 
5pm und 30 um variiert. Der Einfluss auf die Konvertertemperatur 
wird in diesem Abschnitt anhand von Grenz- und Nominalsimulation 
(tg = 15 pm) bei Ip = 1,5 A ermittelt. Die Ergebnisse sind in Abbildung 
5.25 dargestellt: Über den untersuchten Bereich steigt die Temperatur 
zu dickeren Klebeschichten um 7,3 °C. 


(a) ty = 5pm, Tmax = 75,3°C. 


(b) tg = 15 um, Tmax = 79,0°C. 


(c) tg = 30 pm, Tmax = 82,6 °C. 


Abbildung 5.25: Ergebnisse der zeitlich gemittelten Simulation bei Ip = 1,5A mit 
variierender Kleberschichtdicke. Temperaturskala jeweils auf 85°C Maximum. Der 
Temperaturunterschied zwischen tg = 5ym und tg = 30 um beträgt 7,3°C. 
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GEPLANTE LICHTVERTEILUNGEN 


Um zu tiberpriifen, ob die angestrebte Lichtverteilung unterhalb des 
Limits von 150°C bleibt, wurden die in Abbildung 5.26 gezeigten 
Laserleistungsdichten simuliert. 5.26(a) zeigt die zeitlich gemittelte 
Funktion aus der zeitabhangigen Simulation. 5.26(b) zeigt eine Leis- 
tungsdichte, die einen zur Mitte des Scanwegs hin langsamer wer- 
denden MOEMS-Spiegel nachbildet und entspricht einer typischen 
automobilen Fernlichtverteilung mit dem Maximum in der Mitte. Die 
Leistung entspricht in beiden Fallen den maximal zu erwartenden 
12W. 


1045 
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0.5 
0.5 
0 0 
(a) Sechs Laser, wie in der transienten (b) Lichtverteilung mit mittigem Maxi- 
Simulation. mum. 


Abbildung 5.26: Leistungsdichteverteilungen fiir die zeitlich gemittelten Konvertersi- 
mulationen. 
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Man erkennt in den Ergebnissen in Abbildung 5.27, dass die Tem- 
peraturen weit unterhalb der zuvor beschriebenen, kritischen 150 °C 
bleiben. 
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(a) Tmax = 85,2 °C. 
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(b) Tmax = 108,0 °C. 


Abbildung 5.27: Ergebnisse der zeitlich gemittelten Konvertersimulation. Temperatur- 
skala jeweils auf Maximum. Die Leistungsdichteverteilung, die zur Mitte hin zunimmt 
zeigt erwartungsgemäß die höhere Temperatur. Die Temperaturen sind weit entfernt 
vom thermischen Quenching. 
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ZUSAMMENFASSUNG 


Es konnte gezeigt werden, dass die Berechnungsdauer fiir eine zeit- 
lich gemittelte Simulation wesentlich kleiner ist, als die zeitabhan- 
gige Simulation fiir das Erreichen des thermischen Gleichgewichts 
benötigen würde. Außerdem beansprucht sie weniger Arbeitsspei- 
cher. 


Der thermische Widerstand der Konverterunterbaugruppe ist inner- 
halb der simulierten Parameter im Wesentlichen konstant, wie Tabelle 
5.1 zeigt. 


Über die Variation der Kleberschichtdicke konnte gezeigt werden, 
dass mit den gegebenen Randbedingungen innerhalb der Dickenän- 
derung von 25 um eine Temperaturänderung von 7,3°C zu erwarten 


ist. 


Nach der Simulation der angenommenen Pumpleistungsdichtever- 
teilungen aus Abbildung 5.26 sind keine der gezeigten Auswirkun- 
gen des thermischen Quenchings, wie Farbdrift und Lichtstromab- 
nahme, im Einsatz als automobiles adaptives Fernlicht zu erwar- 


ten. 


Weitere Anwendungen für die dreidimensionale keramische Konver- 
tersimulation mittels der Diffusionsnäherung sind LEDs mit kerami- 
schem Leuchtstoff, wenn die Randbedingungen entsprechend ange- 
passt werden oder, wie in [13] argumentiert, die spekularen Rand- 
bedingungen beibehalten werden, da der spekulare Lichteinfall in 
einem hoch streuenden Medium, wie dem keramischen Konverter, 
schnell diffus wird, sodass der Einfluss der Randbedingung als nicht 
signifikant betrachtet werden kann. 


Denkbar ist auch die Diffusionsnäherung zur Simulation thermogra- 
phischer Konverter einzusetzen. Diese dienen zur berührungslosen 
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Temperaturmessung ortsferner Gegenstände. Dazu werden Konverter, 
die meist aus zwei Leuchtstoffen bestehen, auf die zu messende Ober- 
fläche aufgetragen und mittels einer Pumpanregung zum Emittieren 
gebracht [45], [46]. Haben die Leuchtstoffe unterschiedliche thermische 
Quenching-Eigenschaften (siehe Abschnitt 5.3), ändert sich die Farbe 
des emittierten Lichtes mit der Temperatur [45], [46]. Sind die Konver- 
ter hoch streuend, kann man die Diffusionsnäherung anwenden. Man 
muss entsprechend der Leuchtstoffzusammensetzung eine weitere 
„Physik“ für den zweiten Leuchtstoff hinzufügen. 


Um zu überprüfen, ob die Vorhersagen der Simulation mittels der 
Diffusionsnäherung realistische Ergebnisse liefern, werden diese in 
nächsten Kapitel mit Messungen verglichen. 
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Kapitel 6 


KONVERTERMESSUNGEN UND 
VERGLEICH MIT DEM 
DIFFUSIONSSIMULATIONSMODELL 


In diesem Kapitel werden die Simulationsergebnisse aus dem vo- 
rigen Kapitel zunächst mit thermographischen Messungen und im 
Anschluss mit Leuchtdichtemessungen verglichen. 


Die thermischen Messungen werden mit einer Infrarotkamera durch- 
geführt, die Leuchtdichtemessungen mit einer Leuchtdichtekamera. 
Beide Systeme müssen zunächst kalibriert werden. Dazu wird der 
Emissionskoeffizient des Konverters im Infrarotbereich der Thermo- 
graphiekamera bestimmt. Für die Leuchtdichtekamera wird ein Farb- 
vergleich mit einem Spektrometer durchgeführt. Im Anschluss werden 
die jeweiligen Messergebnisse den Simulationsergebnissen gegenüber- 
gestellt. 


6.1 THERMISCHE MESSUNGEN AM LABORAUFBAU 


In diesem Abschnitt soll mittels einer Infrarotkamera die Konverter- 
temperatur unter verschiedenen Randbedingungen gemessen und mit 


dem Simulationsmodell verglichen werden. 
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6.1.1 MESSAUFBAU 


Der Messaufbau ist derselbe wie in Abschnitt 4.3. Abbildung 6.1 zeigt 
eine Fotografie des Aufbaus mit IR-Kamera. Die verwendete Ther- 
mographiekamera ist eine VarioCam HD head der Firma InfraTec 
mit Mikroskopobjektiv. Die Kamera beherbergt einen Bildsensor mit 
1024x768 IR-Pixeln, die mittels MicroScan zu 2048x1536 IR-Pixeln ver- 
rechnet werden können [47]. Weitere Daten sind in Tabelle 6.1 ver- 
zeichnet. 


Tabelle 6.1: Spezifikation der IR-Kamera VarioCam HD head der Firma InfraTec [47]. 


Spektralbereich (7,5 - 14) um 

Detektortyp Ungekühltes Mikrobolometer 
Focal Plane Array 

Detektorformat (1024 x 768) (IR-Pixel) 

Bildformat mit opto- (2048 x 1.536) (IR-Pixel) 


mechanischem MicroScan 

Temperaturmessbereich (-40 - 2000) °C 

Messgenauigkeit +15°C 0< T< 100°C oder 
+15% T < Obzw.T > 100°C 
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Abbildung 6.1: Iteration 10 des allgemeinen Messaufbaus wie in Abbildung 4.7 mit 
IR-Kamera statt Leuchtdichtekamera. 


6.1.2 BESTIMMUNG DES EMISSIONSKOEFFIZIENTEN DES 
KONVERTERS 


Die vergleichende Bestimmung des Emissionskoeffizienten mittels 
einer Infrarotkamera ist in der Literatur gut beschrieben (z. B. [48] 
oder [49]). Der Aufbau ist in Abbildung 6.3 zu sehen und zeigt ein 
Lötbad, in das ein Kupferhalter gesteckt wird, mit dessen Hilfe ei- 
ne Konverterbaugruppe geheizt wird. Eine Hälfte des Konverters 
wurde mit einem Lack mit bekanntem Emissionsgrad lackiert, damit 
ein Temperaturunterschied, wie er beispielsweise durch verschiedene 
thermische Anbindungen hervorgerufen werden könnte, vermieden 
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wird. Im Thermographiebild wurde die Emissivitat im lackierten Teil 
auf den bekannten Wert eingestellt. Nun wurde die Temperatur des 
Lötbads variiert, gewartet bis sich das System im thermischen Gleich- 
gewicht befindet und der Emissionsgrad des unlackierten Teils des 
Konverters im Thermographiebild so angepasst, dass die Temperatu- 
ren von lackiertem und unlackiertem Teil gleich sind. Dies wurde fiir 
Temperaturen von 50 °C bis 250°C in Fiinfgradschritten wiederholt, 
welche nicht immer genau getroffen wurden. Die Ergebnisse sind in 
Abbildung 6.2 dargestellt. Es ergibt sich ein Wert von 0,84 + 0,01 für 
50°C < T < 100°C. 
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Abbildung 6.2: Emissivität des Konverters über seine Temperatur. Im Bereich von 
50°C < T < 100°C ergibt sich ein Wert von 0,84 + 0,01. 


118 


THERMISCHE MESSUNGEN AM LABORAUFBAU 


= 
=, 
ka 

© 
s% | 
= 
= 
= 
= 
= 


Abbildung 6.3: Umfunktioniertes Lötbad zur Erwärmung des Konverters zur Ermitt- 
lung des Emissionskoeffizienten. 


6.1.3 MESSUNGEN MIT SCANNENDEM LASERSPOT 


Die Messung des scannenden Laserspots wurde mit dem in Tabelle 
4.2 charakterisierten 100 %-Spiegel durchgeführt. 


Die maximale Temperatur innerhalb des Messbereichs wurde tiber 
verschiedene Positionen des Konverters entlang der optischen Achse 
und bei verschiedenen Vorwärtsströmen gemessen. 
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Die Ergebnisse sind in Abbildung 6.4 dargestellt. Man erkennt, dass 
die Temperatur über den Verfahrweg nur wenig variiert, wohl aber 
mit dem Vorwärtsstrom. 


x Ip=500mA + Ip=150mA >»  Ir=2000mA 
— DN Ir = 500 mA -DN Ip = 1500 mA --- DN Ir = 2000 mA 


90 


Tmax in °C 


70 F | 
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| | | | 
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Konverterposition entlang der optischen Achse in ym 


Abbildung 6.4: Maximale Temperatur des Konverters über Verschiebung entlang der 
optischen Achse bei verschiedenen Vorwärtsströmen bei Konverterhaltertemperatur 
65 °C im dynamischen Fall mit 100 % Spiegel. 

DN: Simulation mit der Diffusionsnäherung. Für Ir = 1,5 A stimmt die Simulation gut 
mit den Messungen überein. Für Ir = 0,5 A wird die Temperatur in der Simulation 
leicht überschätzt und für Ir = 2,0 A leicht unterschätzt. 


Ebenfalls in Abbildung 6.4 dargestellt sind die Ergebnisse der zeitlich 
gemittelten Simulationen. Für Ir = 1,5 A stimmt die Simulation gut 
mit den Messungen überein. Für Ir = 0,5 A wird die Temperatur in der 
Simulation leicht überschätzt und für Ir = 2,0 A leicht unterschätzt. 
Dies könnte auf zwei Probleme hindeuten: 
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1. Der Konverter zeigt eine stärkere Temperaturabhängigkeit des 
Emissionskoeffizienten, als in Abschnitt 6.1.2 bestimmt. 


2. Die Wärmeleitfähigkeiten der verwendeten Materialien sind so 
stark temperaturabhängig, dass sie nicht näherungsweise als 
konstant hätten angenommen werden dürfen. 


In Abbildung 6.5 sieht man eine Asymmetrie der Temperaturvertei- 
lung. Diese kann durch Schichtdickevariationen im Kleber erklärt 
werden (siehe Abschnitt 5.5.2 und 5.5.4). 


Abbildung 6.5: Thermographieaufnahme des Konverters im dynamischen Fall bei 
2,0 A und Konverterhaltertemperatur 65 °C mit 100 %-Spiegel bei Position 360 pm als 
Beispielbild. Der Messbereich (gelb) R1 gibt den maximalen Temperaturwert aus, der 
blaue Bereich R1 setzt den Emissionskoeffizienten auf 85%. Man erkennt eine leichte 
Asymmetrie in der Temperaturverteilung. 


In Abbildung 6.6 ist ein Schnitt senkrecht zur Scanlinie durch deren 
Mittelpunkt gezeigt. Zu sehen ist die Messung des einfallenden La- 
serstrahls sowie die Temperaturverteilung der Simulation und der 
Messung. Es ist festzustellen, dass im Graphen die Temperatur der 
Messung oberhalb derer der Simulation liegt. Dies deutet auf eine 
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höhere Klebeschichtdicke hin oder, dass der Kleber eine geringere 
thermische Leitfahigkeit hat als angenommen. Dies bedeutet, dass die 
Simulationen aller Wahrscheinlichkeit nach ein wenig zu kühl sind. 
Dies legt den Schluss nahe, dass der Emissionskoeffizient in dem Mess- 
bereich, in dem bei Ir = 0,5 A gemessen wurde, zu niedrig bestimmt 
wurde. 


Außerdem erkennt man in Abbildung 6.6, dass die thermische Ant- 
wort auf den Laser wesentlich breiter ist als die blaue Anregung. Dies 
bedeutet, dass es für ein optimales Kühlsystem wichtig ist, dass der 
Konverter und das Saphirsubstrat wesentlich größer sind als die an- 
regende Pumpquelle, damit die Verlustleistung effizient abgeführt 
werden kann. 


Mit den hier gezeigten Messungen kann abgeschätzt werden, dass mit 
den Temperaturen für die in Abschnitt 5.5.4 simulierten Leistungsdich- 
teverteilungen kein thermisches Quenching im Betrieb als adaptives 
automobiles Fernlicht zu erwarten ist. 
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Abbildung 6.6: Normalisierter Schnitt durch die Mitte des Scanwegs bei Ir = 1,5 A. Die 
thermische Antwort auf den Laser ist wesentlich breiter als die blaue Anregung. Der 
Schnitt der Messung liegt oberhalb des Schnittes der Simulation. 
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6.2 LEUCHTDICHTEMESSUNGEN AM LABORAUFBAU 


In diesem Abschnitt soll mittels einer Leuchtdichtekamera die Leucht- 
dichte unter verschiedenen Randbedingungen gemessen und mit dem 
Simulationsmodell verglichen werden. 


6.2.1 MESSAUFBAU 


Der Messaufbau ist derselbe wie in Abschnitt 4.3. Abbildung 6.7 zeigt 
eine Fotografie des Aufbaus mit Leuchtdichtekamera. Die verwen- 
dete Leuchtdichtekamera ist eine LMK 5 color der Firma Techno- 
Team mit 04314f100 Objektiv. Die Kamera beherbergt einen Bildsensor 
mit 1380 x 1030 Pixeln. Weitere Daten sind in Tabelle 6.2 verzeich- 
net. 


Die Messungen wurden mit einem Neutraldichtefilter durchgeführt, 
um möglichst lange Integrationszeiten und damit die höchst mögliche 
Zahl an Periodendauern in einer Messung zu erhalten. Außerdem 
wurden die Integrationszeiten mit der Scanfrequenz moduliert, um 
Messartefakte auszuschließen. 


Tabelle 6.2: Technische Spezifikation der Leuchtdichtekamera LMK 5 color der Firma 
TechnoTeam [50]. 


Detektortyp CCD 
Detektorformat (1380 x 1030) Pixel 
Integrationszeiten 100 us - 15s 


Kalibriergenauigkeit AL <25% 
Wiederholgenauigkeit AL < 0,1% 
Messgenauigkeit AL < 2,5% für Normlichtart A 
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Abbildung 6.7: Iteration 10 des Messaufbaus zur Bestimmung von Leuchtdichte, Tem- 
peratur und Leistungsdichte. 

1: Leuchtdichtekamera / Thermokamera / Ulbrichtkugel; 2: Leistungsmessgerät; 3: Leis- 
tungsmesskopf (Ophir L50(150)A); 4: Ablenkspiegel; 5: Laserdiode mit Asphare und 
Wasserkiihlung; 6: 2-D Scanspiegel; 7: Konverterhalter / CCD-Kamera 


6.2.2 UNSICHERHEITSBETRACHTUNGEN 


Während der Analyse der Messergebnisse hat sich gezeigt, dass die mit 
der Leuchtdichtekamera gemessenen Farbwerte x und y stark von der 
erwarteten Konversionsgeraden abweichen, vor allem mit den zum 
Schutz der Kamera genutzten Neutraldichtefiltern. Diese Abweichung 
ist mit der in Tabelle 6.2 angedeuteten Kalibrierung mit Normlichtart A 
zu erklären. Dieses Kalibrierlicht entspricht einer Glühfadenlampe mit 
den Farbkoordinaten x = 0,4476 und y = 0,4074 und einer Farbtempe- 
ratur von 2856 K. Normlichtart A weist, im Gegensatz zum Laserscan- 
ner, nur einen sehr kleinen Blauanteil auf. [43] 
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Um diese Abweichung zu bestimmen, wurde mittels eines Spek- 
trometers, an dessen Faserende ein Diffusor montiert wurde, das 
Spektrum des Konverters in einem Meter Abstand aufgenommen, 
um auf der Einheitskugel zu sein. Daraufhin wurden die Farb- 
Leuchtdichteaufnahmen des Konverters angefertigt und die Farbe 
über den gesamten Konverterbereich gemittelt. Die Ergebnisse sind in 
Abbildung 6.8 dargestellt. Die mit der Leuchtdichtekamera gemesse- 
nen Farbwerte weichen von den mit dem Spektrometer ermittelten 
Werte ab. Es ergeben sich abhängig vom benutzen Filter unterschiedli- 
che Werte. 


---ECE-Weißfeld @ Spektrometermessungen x Messungen mit ND3 Filter 
Xx Messungen mit ND5 Filter Linear (Spektrometermessungen) 
v x x 
0.4000 } 
x x 
x xl 
Er 
0.3500 + rn 
coca . 
en 
x = lade 
0.3000 } i se, 
. beers 
0.2500 + e 
0.2000 + 4 + 4 4 
0.2400 0.2600 0.2800 0.3000 0.3200 0.3400 x 


Abbildung 6.8: Farbfehler bei der Leuchtdichtemessung: Die mit der Leuchtdichteka- 
mera gemessenen Farbwerte liegen weit ab der mit dem Spektrometer ermittelten Werte. 
Es ergeben sich abhangig vom benutzen Filter unterschiedliche Werte. 


Die Software der Leuchtdichtekamera lässt es zu, den Farbwert einer 
Aufnahme innerhalb eines bestimmten Messbereichs zu korrigieren. 
Dies wurde fiir den zur Simulation verwendenden Konvertertyp getan. 
Damit stimmt die Farbe fiir diesen einen gemessenen Punkt mit dem 
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(a) Vorher. (b) Nachher. 


Abbildung 6.9: Farbbilder zur Verdeutlichung der Farbfehler in der Leuchtdichtemes- 
sung. Ir = 0,5 A, ND5 Filter. 


Spektrometerwert überein. Alle anderen Farb-Messpunkte weisen wei- 
terhin Fehler auf, wie in Abbildung 6.10 zu sehen ist. Die Fehlerbalken 
stellen lediglich den Abstand zur Spektrometermessung grafisch dar. 
Man erkennt außerdem, dass die Fehler größer werden, je weiter die 
Messungen in Richtung kleiner x- und y-Werte gehen. Ein Vergleich 
der Farbaufnahmen vor und nach der Anpassung der Farbwerte ist in 
Abbildung 6.9 dargestellt. 


Die Berechnung der Leuchtdichte aus der mit dem Spektrometer be- 
stimmten Beleuchtungsstärke E, dem Abstand von Diffusor zu Kon- 
verter r und der Fläche des Leuchtdichtebilds Agg nach Gleichung 
(6.1) kann allenfalls als grobe Abschätzung betrachtet werden, unter 
der Annahme, dass der Konverter Lambert’sch abstrahlt und, dass 
Diffusor und Konverter parallel stehen. Gleichung (6.1) kann nur ange- 
wandt werden, wenn sich der Diffusor auf der Einheitskugel befindet 
und damit r = 1m gegeben ist. 
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Abbildung 6.10: Farbfehler bei der rekalibrierten Leuchtdichtemessung: Die mit der 
Leuchtdichtekamera gemessenen Farbwerte liegen für den einen korrigierten Wert 
(umkreist) nahe am Spektrometerwert. Die restlichen Werte weichen weiterhin von 
den mit dem Spektrometer ermittelten Werte ab. Auch die Konversionsgeraden der 
Leuchtdichtemessungen passen nicht zur Spektrometermessung. 


Vergleicht man die Farbwerte der Simulation aus Abbildung 5.20 mit 
denen der Messung (eingekreist in Abbildung 6.10), so zeigt sich, dass 
die genäherte Farbberechnung der Simulation nicht mit der Messung 
überein stimmt. 


Somit ist eine Aussage über die absolute Farbe und damit auch über 
die absolute Leuchtdichte aus der Simulation nicht machbar. Aussa- 
gen über eine qualitative Leuchtdichte- oder Farbänderung sind aber 
weiterhin möglich. 


Für eine ordentliche Rekalibrierung der Leuchtdichtekamera, wie sie 
beispielsweise in [51] beschrieben ist, fehlte im Projekt die nötige 
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Zeit. Daher sollen im Folgenden lediglich normalisierte Leuchtdichten 
verglichen werden. 


6.2.3 MESSUNGEN MIT SCANNENDEM LASERSPOT 


Die Leuchtdichte mit scannendem Laserspot wurde mit dem in Tabelle 
4.2 beschriebenen 100 %-Spiegel gemessen. Sie wurde an verschiede- 
nen Konverter-Positionen entlang der optischen Achse und bei un- 
terschiedlichen Vorwärtsströmen gemessen. Abbildung 6.11 zeigt ein 
Beispielbild inklusive der Schnittgeraden. 


Die gemessenen, absoluten Maximalwerte der Leuchtdichte unter- 
scheiden sich über den Verfahrweg des Konverters maximal um 5 % 
für Ir = 2,0 A, 6 % für Ir = 1,5 A und 8% für Ir = 0,5 A. Die Messserie 
ist in Tabelle 6.3 dargestellt. 


Abbildung 6.11: Leuchtdichtebild mit Schnittgerade und ohne Skala bei Ir = 1,5 A und 
in Konverterposition 0 pm. 
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Die Abbildungen 6.13 bis 6.15 zeigen Schnitte entlang der Schnitt- 
geraden bei den Vorwärtsströmen Ir = 0,5A, 1,5A und 2,0 A. Man 
erkennt, dass die Breite des Schnitts in 1/e? kaum über den Verfahr- 
weg entlang der optischen Achse variiert. Die gemessenen Werte fiir 
1/e sind in Abbildung 6.12 dargestellt. Die Variation ist gering, war 
aber nach den Ergebnissen aus Abschnitt 4.3 kleiner erwartet wor- 
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Abbildung 6.12: Scanlinienbreitenauswertung in 1/e für verschiedene Vorwärtsströme 
über der Konverterposition entlang der optischen Achse. 
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Tabelle 6.3: Messserien zur Bestimmung der maximalen Leuchtdichteabweichung in- 
nerhalb einer Vorwärtsstrom-Messreihe. 


Position pm E rel (0.5 A) I riesg rel (1.5 A) neg rel (2.0 A) 


-450 94% 94% 98% 
-405 92 % 96 % 96 % 
-360 94% 95% 95% 
-315 100 % 97% 98% 
-270 97% 97% 97% 
-225 99% 99% 99 % 
-180 100 % 100 % 98 % 
-135 97 % 98 % 98 % 
-90 98% 100 % 100 % 
-45 99 % 99 % 99 % 
0 100 % 100 % 99 % 
45 96 % 97 % 98 % 
90 97 % 98 % 98 % 
135 99 % 98 % 98 % 
180 99 % 99 % 99 % 
225 99 % 99% 99% 
270 98% 99% 98% 
315 100 % 98% 98% 
360 99% 97% 98% 
405 95% 97% 97% 
450 94% 96 % 97 % 
max. Abweichung 8% 6% 5% 
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Abbildung 6.13: Normalisierter Leuchtdichteschnitt durch die Mitte des Scanwegs 
bei Ir = 0,5 A. Die Scanlinienbreite variiert maximal um 66 um über den Verfahrweg 
entlang der optischen Achse. Die Einträge in der Legende tragen die Einheit pm. 
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Abbildung 6.14: Normalisierter Leuchtdichteschnitt durch die Mitte des Scanwegs 
bei Ir = 1,5 A. Die Scanlinienbreite variiert maximal um 67 um über den Verfahrweg 
entlang der optischen Achse. Die Einträge in der Legende tragen die Einheit pm. 
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Abbildung 6.15: Normalisierter Leuchtdichteschnitt durch die Mitte des Scanwegs 
bei Ir = 2,0 A. Die Scanlinienbreite variiert maximal um 90 um über den Verfahrweg 
entlang der optischen Achse. Die Einträge in der Legende tragen die Einheit pm. 
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Abbildung 6.16: Normalisierter Leuchtdichteschnitt durch die Mitte des Scanwegs bei 
verschiedenen Vorwärtsströmen in Konverterposition 0 pm. Die Scanlinienbreite variiert 


maximal um 38 pm. 
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6.2.4 VERGLEICH DER SCANLINIENBREITE VON LEUCHTDICHTE 
UND BLAUER PUMPANREGUNG 


Vergleicht man die Breite der Scanlinie des Lasers mit der der Leucht- 
dichte, wie in Tabelle 6.4 geschehen, so fallt auf, dass die Leuchtdichte 
mit einem Faktor von 3,7-4,0 breiter in 1/e ist. Dies limitiert den 
möglichen Kontrast zwischen zwei Scanlinien und während des Aus- 
blendens des Gegenverkehrs. 


Tabelle 6.4: Vergleich der Mittelwerte der gemessenen 1/e?-Breiten des Lasers und der 
Leuchtdichte und der sich daraus ergebende Vergrößerungsfaktor. 


Vorwärtsstrom 1/e Laser 1/e? Leuchtdichte Verbreiterungs- 


in A in pm in ym faktor 
0,5 169 671 3,7 
15 173 681 3,9 
2,0 172 687 4,0 


Ebenso konnte mit den Messungen aus 4.3 erwartet werden, dass 
die Scanlinienbreite der Leuchtdichte wenig mit dem Vorwärtsstrom 
variiert. Abbildungen 6.12 und 6.16 zeigen den geringen Einfluss des 
Vorwartsstroms auf die Scanlinienbreite. 


Da die Scanlinienbreite weder in der Messung noch in der Simulation 
(siehe Abschnitt 5.5.4) wesentlich mit dem Vorwärtsstrom oder der 
Konverterposition variiert, genügt es ein Schnittbild zum Vergleich 
von Simulation und Messung heran zu ziehen: In Abbildung 6.17 ist 
zu erkennen, dass der Verlauf der Leuchtdichte von Simulation und 
Messung oberhalb von 1/e? sehr gut übereinstimmt. Unterhalb von 1/e? 
liegt die Messung leicht oberhalb der Simulation. Das bedeutet eine 
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weitere Verschlechterung des Kontrastes zwischen zwei Scanlinien 


und während des Ausblendens des Gegenverkehrs. 


Um zu beurteilen, ob dieser Effekt an Mehrfachreflexionen im Umlenk- 
spiegel liegt (Position (4) in Abbildung 6.7), wurde eine Simulation mit 
dem zweiten Reflex aus der Messung mit dem 90/10 Strahlteiler durch- 
geführt. Das Ergebnis ist in Abbildung 6.18 gezeigt. Man erkennt, dass 
ein zweiter Reflex zu nah am primären Laserstrahl verortet wäre und 
sich schon oberhalb von 1/e? eine deutliche Abweichung zeigen müss- 
te. Damit kann eine Mehrfachreflexion im Umlenkspiegel als Grund 
für die Abweichung ausgeschlossen werden. 


Als weitere, sehr wahrscheinliche Fehlerquelle, kommen die in der 
Simulation mit der Diffusionsnäherung getroffenen Annahmen in 
Frage. Durch die Scanbewegung trifft der Laserstrahl nicht senk- 
recht auf den Saphir. Damit kommt es im Saphir zu Mehrfachrefle- 
xionen des Laserstrahls. Diese können die Konverterantwort aufwei- 
ten. 
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Abbildung 6.17: Normalisierter Schnitt durch die Mitte des Scanwegs bei Ir = 1,5 A. 
Die Scanlinienbreite der Messung variiert im Vergleich zur Simulation oberhalb von 1/e? 
nur wenig. Unterhalb von 1/e? liegt die Messung über der Simulation. 
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Abbildung 6.18: Normalisierter Schnitt durch die Mitte des Scanwegs bei Ir = 1,5A 
inklusive des zweiten Reflexes aus den Messungen aus Abschnitt 4.3. 
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6.3 ZUSAMMENFASSUNG 


Für die thermischen Messungen in Abschnitt 6.1 wurde zunächst der 
Emissionskoeffizient des Konverters bestimmt. In der Folge konnten 
die Simulationsergebnisse mit Thermografiemessungen verglichen 
werden. Es zeigte sich, dass in den Simulationen der Kleber wahr- 
scheinlich ein wenig zu dünn oder mit zu großer Warmeleitfahigkeit 
angenommen wurde, sodass die Simulation bei Ir = 1,5 A maximal 
2°C und bei Ir = 2,0 A höchstens 5 °C kühler ist, als die Messung. Die 
Ergebnisse bei Ir = 0,5 A legen nahe, dass der Emissionskoeffizient 
des Konverters für niedrige Temperaturen als zu hoch bestimmt wur- 
de. Nach den thermographischen Messungen ist davon auszugehen, 
dass in der Anwendung als automobiles, adaptives Fernlicht kein 
thermisches Quenching und damit ein sicherer Betrieb zu erwarten 
ist. 


Für die Leuchtdichtemessungen in Abschnitt 6.2 wurde zunächst eine 
grobe Farbkalibrierung durchgeführt. Vergleicht man die Farbwerte 
der Simulation aus Abbildung 5.20 mit denen der Messung (einge- 
kreist in Abbildung 6.9), so zeigt sich, dass die genäherte Farbberech- 
nung der Simulation nicht mit der Messung überein stimmt. Somit ist 
eine Aussage über die absolute Farbe und damit auch über die abso- 
lute Leuchtdichte aus der Simulation nicht machbar. Aussagen über 
eine qualitative Leuchtdichte- oder Farbänderung sind aber weiterhin 
möglich. 


Vergleicht man die Breite der Scanlinie des Lasers mit der der Leucht- 
dichte, wie in Tabelle 6.4 geschehen, so fällt auf, dass die Leuchtdichte 
mit einem Faktor von 3,7-4,0 breiter in 1/e ist. Dies limitiert den 
möglichen Kontrast zwischen zwei Scanlinien und während des Aus- 
blendens des Gegenverkehrs. 


138 


ZUSAMMENFASSUNG 


In Abbildung 6.17 ist zu erkennen, dass der Verlauf der Leuchtdich- 
te von Simulation und Messung oberhalb von 1/e? sehr gut überein- 
stimmt. Unterhalb von 1/e liegt die Messung leicht oberhalb der Si- 
mulation. Das bedeutet eine weitere Verschlechterung des Kontrastes 
zwischen zwei Scanlinien und während des Ausblendens des Gegen- 
verkehrs. Eine Mehrfachreflexion im Umlenkspiegel konnte als Grund 
für die Abweichung ausgeschlossen werden. Als sehr wahrscheinliche 
Fehlerquelle kommen die in der Simulation mit der Diffusionsnähe- 
rung getroffenen Annahmen in Frage: Durch die Scanbewegung trifft 
der Laserstrahl nicht senkrecht auf den Saphir. Damit kommt es im 
Saphir zu Mehrfachreflexionen des Laserstrahls. Diese können die 


Konverterantwort aufweiten. 
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ZUSAMMENFASSUNG UND AUSBLICK 


7.1 ZUSAMMENFASSUNG 


Diese Arbeit hat in Kapitel 2 gezeigt, dass die thermischen Wider- 
stande der zur Verfiigung stehenden Laserdioden mit der Temperatur, 
jedoch nicht wesentlich untereinander variieren. Lediglich die Diode 
im TO56-Gehäuse sticht heraus, da sie bei Erzeugung von 1 W opti- 
scher Leistung 3,5 W thermische Verlustleistung erzeugt und damit 
ineffizienter arbeitet als die restlichen Varianten. Die drei anderen 
Laserdioden generieren bei 2W Ausgangsleistung nur 5,7W Abwär- 
me. Betrachtet man sie einzeln, ohne Halter und Anbindung an ein 
Kühlsystem, könnte man aus thermischer Sicht jede der Varianten 


wählen. 


Je höher die mögliche maximale Interface-Temperatur zwischen La- 
sermodul und dessen Kühlsystem ist, desto kleiner und günstiger 
kann dieses ausfallen. In Abschnitt 3.1 wurde gezeigt, dass die Lösung 
mit SMD-Laserdiode und TO90-Laserdiode gleichwertig sind und die 
Lösungen mit MDP und TO56 Lasern dagegen abfallen. Wegen der 
besseren Verfügbarkeit wurde die TO90-Laserdiode ausgewählt. Die 
thermografische Messung des Prototypen des Laserdiodenmoduls 
passt sehr gut zu den Simulationen, wie in Abschnitt 3.2 dargelegt 


wurde. 


Das Interface des Lasermoduls zum Kühlsystem ist wie folgt beschrie- 
ben worden: Jede TO90 Laserdiode trägt eine thermische Last von 
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4,3W. Am Interface zwischen Kühlkörper und Lasermodul dürfen bei 
dieser Last und Lastverteilung maximal 79 °C anliegen. Die zur Wär- 
meübertragung nutzbare Fläche beträgt 34,2 cm?. 


Thermomechanische Simulationen der TO90 Laserdiode in Abschnitt 
4.1 und der TO90-Asphären-Unterbaugruppe in Abschnitt 4.2 zeigen 
auf, dass sich die Facettenposition mit der Temperatur ändert und, 
dass die temperaturbedingte Ausdehnung auf den Linsenhalter ver- 
gleichsweise geringe Auswirkungen hat. Axiale Schnitte in der fast 
und slow axis des Laserprofils veranschaulichen den Einfluss der Tem- 
peratur detaillierter. Die effektive Verschiebung des Laserprofils in der 
fast axis beträgt 84 um und in der slow axis 98 um. Die Defokussierung 
in Richtung der asphärischen Linse führt zu einer ungewünschten Ver- 
ringerung der Leistungsdichte. 


Der Einfluss des Vorwärtsstroms auf die Fokusposition des hier ver- 
wendeten optischen Systems wurde in Abschnitt 4.3 gemessen. Es 
stellte sich heraus, dass die Fokusposition über den Vorwärtsstrom 
um mehrere Millimeter variiert. Für Anwendungen in denen höchste 
Leuchtdichten gefordert sind, muss die Asphärenposition genau auf 
den Vorwärtsstrom angepasst sein. 


Eine solche Nutzung ist beispielsweise ein LARP-Zusatzfernlicht, das 
Sichtweiten von über 600 m ermöglicht. Weitere LARP-Anwen-dungen 
sind Konverterräder, wie sie in Projektoren oder zur Endoskopie ein- 
gesetzt werden, bei denen Leuchtdichten von über 500 cd - mm? oder 
gar über 1000 cd - mm? gefordert sind, was selbst die leistungsstarks- 
ten LEDs heute nicht zu leisten vermögen [33]. 


Betrachtet man zusätzlich die Änderung der Facettenposition über 
die Temperatur, müssen in der Montage eines jeden einzelnen TO90- 
Submounts die Laserdiode und die Asphäre unter Betriebsbedingun- 
gen, das bedeutet unter den im Fahrzeug vorherrschenden Bedingun- 
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gen, bei Betriebsstrom (hier: Ir = 2A) und bei Betriebstemperatur 
(hier: Te = 70°C) gegeneinander justiert werden. 


Diese Erkenntnis führte dazu, dass für die Nachfolgeprojekte der p- 
LARP-Scheinwerfer eine aktive Justage von Asphäre gegen Facette in 
der Fertigung etabliert wurde. Dadurch konnten die Leuchtdichte er- 
höht und die Leuchtdichte- und Spotpositionstoleranz zweier u-LARP- 
Gen2-Module verringert werden. Dadurch lassen sich kleinere Optiken 
und/oder größere „Lichtreichweiten” erzielen. 


In Abschnitt 4.3 wurde auch die Auswirkung des Vorwärtsstroms 
auf die Scanlinienbreite und die Leistungsdichte bei sich entlang der 
optischen Achse bewegendem Detektor gemessen. Es stellte sich her- 
aus, dass sich die Leistungsdichte innerhalb des Scanwegs kaum über 
die Verschiebung entlang der optischen Achse ändert. Darin könn- 
te ein weiterer Vorteil dieser Anordnung liegen, denn man könnte 
die punktuelle Leistungsdichte über eine Defokussierung reduzie- 
ren, um optisches Quenching des Leuchtstoffs zu verringern (siehe 
Abschnitt 5.3). Dies müsste in weiteren Messungen untersucht wer- 
den. 


Eine weitere Folge der sich kaum ändernden Leistungsdichte und 
Scanlinienbreite über die Verschiebung entlang der optischen Achse ist, 
dass in der Simulation nur eine geringe Änderung der Temperatur für 
die gemessenen Scanlinien innerhalb einer Vorwärtsstrom-Messreihe 


zu erwarten ist. 


Möchte man ein optisches Konzept mit möglichst kleinen, überein- 
ander angeordneten Scanlinien realisieren, muss man die slow axis 
sehr genau entlang der Scanlinie ausrichten. Selbst bei kleinen Ab- 
weichungen dominiert die Breite der slow axis über die Genauigkeit 
der Konverterposition entlang der optischen Achse. Setzt man im 
Gegensatz dazu die fast axis mit der Scanrichtung gleich, ergibt sich 
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ein gegen Rotation des Laserspots wesentlich unempfindlicheres Sys- 


tem. 


In Kapitel 5 wurde zum ersten Mal die Anwendung der Diffusionsnä- 
herung der Strahlungstransportgleichung auf einen dreidimensiona- 
len, nicht rotationssymmetrischen, keramischen Leuchtstoff mit über 
der Zeit örtlich veränderlicher Laseranregung gezeigt. Dazu wurden 
zunächst die Grundlagen des thermischen und optischen Quenchings 
sowie der Diffusionsnäherung erläutert und deren Einbindung in Com- 
sol Multiphysics® beschrieben. Es hat sich herausgestellt, dass eine 
unphysikalische Dämpfung einzuführen ist, um sinnvolle Ergebnisse 
zu erhalten. Die Simulationszeiten für die transiente Simulationen sind 
jedoch so hoch, dass sie mit den zur Zeit zur Verfügung stehenden 
Rechnern nicht wirtschaftlich zu berechnen sind. 


In Abschnitt 5.5.4 wurden zeitlich gemittelte Simulationen mit den 
zuvor gemessenen Leistungsdichteverteilungen durchgeführt. Diese 
Berechnungen liefen schneller ab als die zeitabhängigen Simulationen. 
Ferner stellte sich heraus, dass der thermische Widerstand nicht von 
der Leistung abhängt. Außerdem wurde eine parametrische Studie mit 
veränderlicher Konverterkühlkörpertemperatur durchgeführt, um die 
Auswirkungen des thermischen Quenchings zu verdeutlichen. Die Er- 
gebnisse wurden mit thermografischen und Leuchtdichtemessungen 
verglichen. 


Nach der Simulation der im Scheinwerfer erwarteten Pumpleistungs- 
dichteverteilungen sind keine der gezeigten Auswirkungen des ther- 
mischen Quenchings wie Farbdrift und Lichtstromabnahme im Einsatz 
als automobiles, adaptives Fernlicht zu erwarten. 


Für die thermischen Messungen in Abschnitt 6.1 wurde zunächst der 
Emissionskoeffizient des Konverters bestimmt. In der Folge konnten 
die Simulationsergebnisse mit Thermografiemessungen verglichen 
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werden. Es zeigte sich, dass in den Simulationen der Kleber wahr- 
scheinlich ein wenig zu dünn oder mit zu großer Warmeleitfahigkeit 
angenommen wurde, sodass die Simulation bei Ir = 1,5 A maximal 
2°C und bei Ir = 2,0 A höchstens 5 °C kühler ist, als die Messung. Die 
Ergebnisse bei Ir = 0,5 A legen nahe, dass der Emissionskoeffizient 
des Konverters für niedrige Temperaturen als zu hoch bestimmt wurde. 
Auch nach den thermografischen Messungen ist davon auszugehen, 
dass in der Anwendung als automobiles, adaptives Fernlicht kein 
thermisches Quenching und damit ein sicherer Betrieb zu erwarten 
ist. 


Für die Leuchtdichtemessungen in Abschnitt 6.2 wurde zunächst ei- 
ne Farbkalibrierung durchgeführt. Vergleicht man die Farbwerte der 
Simulation mit denen der Messung, so zeigt sich, dass die genäherte 
Farbberechnung der Simulation nicht mit der Messung übereinstimmt. 
Somit ist eine Aussage über die absolute Farbe und damit auch über 
die absolute Leuchtdichte aus der Simulation nicht machbar. Aussa- 
gen über eine qualitative Leuchtdichte- oder Farbänderung sind aber 
weiterhin möglich. 


Vergleicht man die Breite der Scanlinie des Lasers mit der der Leucht- 
dichte, so fällt auf, dass die Leuchtdichte mit einem Faktor von 
3,7-4,0 breiter in 1/e? ist. Dies limitiert den möglichen Kontrast 
zwischen zwei Scanlinien und während des Ausblendens des Gegen- 
verkehrs. 


Der Verlauf der Leuchtdichte von Simulation und Messung stimmen 
oberhalb von 1/e sehr gut überein. Unterhalb von 1/e liegt die Mes- 
sung leicht oberhalb der Simulation. Das bedeutet eine weitere Ver- 
schlechterung des Kontrastes zwischen zwei Scanlinien und während 
des Ausblendens des Gegenverkehrs. Eine Mehrfachreflexion im Um- 
lenkspiegel konnte als Grund für die Abweichung ausgeschlossen 
werden. Als sehr wahrscheinliche Fehlerquelle kommen die in der 
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Simulation mit der Diffusionsnäherung getroffenen Annahmen in Fra- 
ge. Durch die Scanbewegung trifft der Laserstrahl nicht senkrecht 
auf den Saphir. Damit kommt es im Saphir zu Mehrfachreflexionen 
des Laserstrahls. Diese können die Konverterantwort aufweiten. Die- 
se Mehrfachreflexionen könnten zukünftig mittels eines ray tracing 


® 


Programms simuliert und in Comsol Multiphysics” eingelesen wer- 


den. 


7.2 AUSBLICK 


Ein besonders interessanter Aspekt dieser Arbeit war die erstmalige 
transiente Simulation eines keramischen Leuchtstoffs mit scannender 
Pumplichtverteilung mittels der Diffusionsnäherung. Problematisch 
war die lange Berechnungsdauer von vier Wochen für zehn Perioden- 
dauern. 


Während die Dichte von Transistoren auf Silizium nach dem Moo- 
re’schen Gesetz weiter zunimmt (wenn auch langsamer als in der 
Vergangenheit), haben Leistungsdichtebeschränkungen die Taktrate 
von Prozessoren seit mehr als einem Jahrzehnt begrenzt [52]. Dies 
hat effektiv jeden Fortschritt in der Verarbeitungskapazität eines ein- 
zelnen CPU-Kerns aufgehalten [52]. Stattdessen ergibt sich eine Er- 
höhung der potenziellen Verarbeitungskapazität aus einer Erhöhung 
der Kernanzahl von CPUs und breiten CPU-Registern [52]. Alternati- 
ve Verarbeitungsarchitekturen sind üblicher geworden. Diese reichen 
von der Mehrkernarchitektur, die auf Standard-x86_64-Kernen basiert, 
bis zu zahlreichen Alternativen wie GPUs. Bei GPUs kann ein viel 
größerer Teil des Chips für arithmetische Berechnungen verwendet 
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werden. Jedoch gehen damit Schwierigkeiten bei der Programmie- 
rung und der Speicherhandhabung einher [52]. Mit diesen Aussichten 
könnte die Berechnungsdauer für die transiente Simulation reduziert 


werden. 


Zusätzlich zu einem scannenden Lasersystem lässt sich die zeitab- 
hängige Simulation mittels der Diffusionsnäherung auch auf Phos- 
phorräder, wie sie in Projektoren vorkommen oder in der Endoskopie 
Verwendung finden, anwenden. Dazu müssen die entsprechenden 
Streukoeffizienten und Reflexivitäten gemessen werden, sollte es sich 
um in Silikon oder Wasserglas eingebettete Konverter handeln. Diese 
müssen darauf folgend auf ihre Tauglichkeit für die Diffusionsnä- 
herung untersucht und entsprechend ins Modell eingesetzt werden. 
Hier lassen sich eventuell Symmetrieeigenschaften nutzen, um die 


Berechnungszeit zu verkürzen. 


Ein weiteres Anwendungsgebiet könnte die Analyse von keramischen 
Konvertern auf LEDs im pulsweitenmodulierten Betrieb sein. Da hier 
sehr ähnliche Konverter Verwendung finden, ist lediglich die Rand- 
bedingung der spekularen Pumpanregung zu ändern. In [13] wird 
dagegen argumentiert, dass die spekularen Randbedingungen beibe- 
halten werden können, da der spekulare Lichteinfall in einem hoch 
streuenden Medium, wie dem keramischen Konverter, schnell diffus 
wird, sodass der Einfluss der Randbedingung als nicht signifikant 
betrachtet werden kann. Auch bei der Simulation von LEDs lassen 
sich die Symmetrieeigenschaften der LED nutzen, um Simulationszeit 


zu sparen. 


Betrachtet man auch Alterungsmechanismen des Konverters und passt 
die Streukoeffizienten in der Simulation entsprechend an, kann man 
auch gealterte, hoch streuende, keramische Konverter, egal in welcher 
Anwendung, simulieren, und die Unterschiede zwischen Tag eins und 
am Ende der Lebensdauer herausarbeiten. 
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Denkbar ist auch die Diffusionsnäherung zur Simulation thermografi- 
scher Konverter einzusetzen. Diese dienen zur berührungslosen Tem- 
peraturmessung ortsferner Gegenstände. Dazu werden Konverter, die 
meist aus zwei Leuchtstoffen bestehen, auf die zu messende Ober- 
fläche aufgetragen und mittels einer Pumpanregung zum Emittieren 
gebracht [45], [46]. Haben die Leuchtstoffe unterschiedliche thermische 
Quenching-Eigenschaften (siehe Abschnitt 5.3), ändert sich die Farbe 
des emittierten Lichtes mit der Temperatur [45], [46]. Sind die Konver- 
ter hoch streuend, kann man die Diffusionsnäherung anwenden. Man 
muss entsprechend der Leuchtstoffzusammensetzung eine weitere 
„Physik“ für den zweiten Leuchtstoff hinzufügen. 


In [33] wird ein Modell diskutiert, das das optische Quenching 
theoretisch beschreibt. Dieses Modell ließe sich in die bestehende 
Simulation der Diffusionsnäherung einbauen. Damit ließen sich 
genauere Vorhersagen über die zu erwartenden Leuchtdichteän- 
derungen bei verschiedenen Leistungsdichten der Pumpanregung 


machen. 


Seit Mitte 2016 bietet Comsol Multiphysics® ein ray tracing add on 
an. Dies behebt ein Problem, das es zuvor zeit- und arbeitsaufwendig 
machte, ray tracing mit thermischen Simulationen zu verbinden: Man 
musste zunächst eine ray tracing Simulation in Programm A laufen 
lassen, und sich die Verluste ausgeben lassen. Diese fügte man in Pro- 
gramm B ein und ließ eine Thermosimulation laufen. Die Ergebnisse 
fügte man wieder in Programm A ein, um die Auswirkungen auf die 
optische Simulation zu analysieren. Danach mussten die Ergebnisse 
wieder in Programm B eingefügt werden, und so weiter. Dies vermei- 
det die Diffusionsnäherung, da hier der optische und thermische Teil 
selbstkonsistent gelöst werden. 


Mittels des neuen Comsol Multiphysics® Moduls könnte man den 
optischen Teil der Simulation lösen und gleichzeitig und gekoppelt 
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den thermischen Teil in einer FEM-Simulation. Es gibt bereits erste 
Arbeiten, die die Kombination aus ray tracing und thermischer FEM- 
Simulation nutzen. Beispielsweise in der Architektur [53] oder zur 
Betrachtung von thermischen Solarkollektoren [54]. In [55] wird die 
Nutzung zur Berechnung des thermal lensing Effekts in Laserdioden 
beschrieben, in der neben dem ray tracing und dem Wärmeleitungs- 
Modul auch das Modul zur mechanischen Simulation benutzt wird. 
In einer Konvertersimulation liefsen sich damit auch Vorhersagen tiber 
die absolute Leuchtdichte machen, da keine Näherung der Strahlungs- 
transportgleichung gemacht werden miisste. 
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Anhang A 


GRAPHEN AUS DEN VORLAUFIGEN 
LASER-DATENBLATTERN 


Auf den folgenden Seiten sind die Graphen aus den vorlaufigen Daten- 
blättern der untersuchten Laserdioden dargestellt. 


169 


ANHANG 


A.1 TO56 LASERDIODE 
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(a) Temperaturabhängige optische Leis- (b) Vorwartsspannung über Vorwärts- 


tung über Vorwärtsstrom [11]. strom [11]. 
| | | | 
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& 12 | Kos 
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(c) Optische Maximalleistung über Gehäu- (d) Temperaturabhängiger Schwellstrom 
setemperatur [11]. [11]. 


Abbildung A.1: Graphen aus dem vorläufigen Datenblatt der TO56 Laserdiode zur 
Berechnung der thermischen Leistung. 
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A.2 TO90 LASERDIODE 


— — 40°C —— 25°C —— 85°C — —40°C —— 25°C —— 85°C 
— 0°C — 55°C — 105°C —— 0°C — 55°C — 105°C 
T T T T T T T T 
3r 5 - 
= > 
a2 £ 
& a i 
a if 
0 LL | ji | | 3 | | i | | 
0 1 2 3 0 1 2 3 
Ir in A Ig inA 


(a) Temperaturabhängige optische Leis- 
tung über Vorwärtsstrom [12]. 
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(c) Temperaturabhängige optische max. 


Leistung über Gehäusetemperatur 


[12]. 


[12]. 


(d) Temperaturabhangiger Schwellstrom 


Abbildung A.2: Graphen aus dem vorläufigen Datenblatt der TO90 Laserdiode zur 
Berechnung der thermischen Leistung. 


Anhang B 


MATERIALDATEN AUS DER 
SOLIDWORKS DATENBANK 


Auf den folgenden Seiten sind die Materialdaten aus der SolidWorks 
Datenbank dargestellt. 
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B.1 KUPFER 


Path: Solids Pre-Defined \ Metals 
Density: 8960 kg - m7? 
Conductivity type: Isotropic 
Electrical conductivity: Conductor 
Radiation properties: No 
Melting temperature: Yes 
Temperature: 1083,05 °C 
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aso ¥ 
= 400 E 25000 
» 350 z 20000 
S 300 £ 
5 250 & 15000 
g 20 g 
S 150 3 10000 
# 100 È 
74 5 5000 
50 E 
0 0 
500 0 500 1000 1500 -500 0 500 1000 1500 
Temperature["C] Temperature[*C] 
(a) Spezifische Wärmekapazität über (b) Wärmeleitfähigkeit über Temperatur. 
Temperatur. 
1.2E-07 
1.0E-07 
E 
E 8.0E-08 
pa 
[eJ 
= 6.0E-08 
E 
2 406-08 
č 
2.0E-08 
0.0E+00 
-500 0 500 1000 1500 


Temperature[*C] 


(c) Spezifischer elektrischer Widerstand 
über Temperatur. 


Abbildung B.1: Graphen aus der Materialdatenbank von SolidWorks zum Material 
Kupfer. 
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B.2 ALUMINIUM 


Path: Solids Pre-Defined \ Metals 
Density: 2688,90 kg - m7’ 
Conductivity type: Isotropic 
Electrical conductivity: Conductor 
Radiation properties: No 
Melting temperature: Yes 
Temperature: 660,25 °C 
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„_ 1200 = 30000 
* 1000 > 
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Temperature[’C] 
(c) Spezifischer elektrischer Widerstand 
über Temperatur. 


Abbildung B.2: Graphen aus der Materialdatenbank von SolidWorks zum Material 
Aluminium. 
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B.3 STAHL 


Path: Solids User Defined\3297700swm000.SLDASM \ 
2,0W 75% Andere Ubergange 

Density: 7700kg : m~’? 

Specific heat: 460J - kg7! - K7! 

Conductivity type: Isotropic 


Thermal conductivity: 30W -m71 . K7! 
Electrical conductivity: Conductor 


8.0E-07 
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Temperature[°C] 


Abbildung B.3: Spezifischer elektrischer Widerstand tiber Temperatur aus der Material- 
datenbank von SolidWorks zum Material Stahl. 
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B.4 WARMEAUSDEHNUNGSKOEFFIZIENTEN 


Tabelle B.1: Materialdaten bei 25 °C fiir die CTE-Simulation der TO90 Laserdiode [13]. 


Schicht Material Wärmeausdehnungskoeff. 


a/c PA 
Gehäuse Stahl 11 
Balkon Cu 24 
Submount 6H-SiC 4,2 / 4,68 
Die GaN 5,59 / 3,17 
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